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ABSTRACT 

We study systematically which features in the cosmic microwave background (CMB) and large-scale structure 
(LSS) probe various inhomogeneous properties of the dark sectors (including neutrinos, dark matter, and dark 
energy). We stress, and quantify by simple formulas, that the primary CMB anisotropics are very susceptible to 
the gravitational potentials during horizon entry, less at recombination. The CMB thus allows us to scan <fr + 
and the underlying dark kinetics for all redshifts z ~ 1 - 10 5 . LSS, on the other hand, responds strongest to 4> 
at low redshifts. Dark perturbations are often parameterized by the anisotropic stress and effective sound speed 
(stiffness). We find that the dark anisotropic stress and stiffness influence the visible species at the correspondingly 
early and late stages of horizon entry, and affect stronger respectively the CMB and LSS. The CMB yet remains 
essential to probing the stiff perturbations of light neutrinos and dark energy, detectable only during horizon entry. 
The clustering of dark species and large propagation speed of their inhomogeneities also map to distinctive features 
in the CMB and LSS. Any parameterization of the signatures of dark kinetics that assumes general relativity can 
effectively accommodate any modified gravity (MG) that retains the equivalence principle for the visible sectors. 
This implies that formally the nonstandard structure growth or $ /\]/ ratio, while indicative, are not definitive MG 
signatures. The definitive signatures of MG may include the strong dependence of the apparent dark dynamics on 
visible species, its superluminality, and the nonstandard phenomenology of gravitational waves. 

Subject headings: cosmology: theory — cosmic microwave background — large-scale structure of universe — 
dark energy — modified gravity 



1. INTRODUCTION 

Galactic and cluster dynamics, cosmic structure, type la su- 
pernovae, the cosmic microwave background (CMB), and the 
primordial abundances of light elements provide solid evidence 
that dark sectors constitute a significant energy fraction of the 
universe at any accessible redshift z < 10 1() . At all correspond- 
ing cosmological epochs the nature of abundant dark species, 
coupled to photons and baryons only by gravitation, is partly or 
entirely uncertain. 

The mainstream analyses of cosmological data usually as- 
sume the minimal neutrino sector, non-interacting cold dark 
matter (CDM), and dark energy represented by a canoni- 
cal scalar field (quintessence). These assumptions are rea- 
sonable for interpreting the available data, yet none of them 
can be taken for granted. For example, new light weakly 
interacting particles commonly appear in high-energy mod- 
els. In some models, even the standard neu trinos recou- 
ple to each other or to addition a l light field s dChacko et all 
2004, 120051; iBeacom et all 12004: Okuil 120051: iGrossman et ail 



2005) at redshifts at which the decoupled component of ra- 



diati on gravitationally a ffects the CMB and cosmic structure 
dHu & Sugivamalll996t iBashinskv & Seliakl l2004t fHannestadl 
2005: iBell et alj|2006l) . Various alternatives to cold dark mat- 
ter hav e been suggested as we ll . These include warm dark 
matter dBlumenthal et al.l Il982t lOlive & Turner! 1 19821) self- 
interacting dark matte r ( Carlson et al J 1 992 : deLaix et al J 1 9951: 
Spergel & Steinhardt 2000), or modified gravity dMilgroml 
1983pBekensteinl2004l:ISkordis et alj2006tlDodelson & Liguoril 
2006). The viability of such scenarios remains an intrigu- 
ing question. Quintessence models are convenient for quan- 
titatively constraining dark energy parameters by data. Yet 
quintessence is not readily motivated by particle physics, where 
it is difficult to naturally achieve the required shallowness of the 
field potential. On the other hand, many alternatives have been 



proposed whose inhomogeneous kinetics, and hence cosmo- 
logical signatures, cannot be mimicked by quintessence with 
any background equation of state w(z) . (For a comprehensi ve 
review of dark energy models see, e.g., |Copeland et al]|2006l) . 

Fortunately, cosmological observations themselves can test 
these assumptions by revealing not only the dark species' mean 
density and pressure but also the kinetics of their inhomo- 
geneities. The goal of this paper is to map various inhomo- 
geneous kinetic properties of the dark sectors (or deviations 
from Einstein gravity) to the observable characteristics of CMB 
and cosmic structure. Dark species influence the visible mat- 
ter by affecting both the background expansion and metric per- 
turbations. Of the two mechanisms, the perturbations, albeit 
demanding better statistics for useful constraints, encode many 
more independent clues about the dark universe by offering new 
information at every spatial scale k. The following three ex- 
amples show the importance of this information, absent in the 
background equation of state w(z). 

1.1. Examples of the value of dark perturbations 
Nature of dark energy 

The first example is the most challenging problem in today's 
cosmology — the nature of dark energy. The constraints on the 
dark energy background equation of state w = /?de/Aie are tight- 
ening around the value -1, consistent with a cosmological con- 
stant. Analyses that combine the current data from the CMB, 
large scale structure (LSS), Lyman-a forest, and supernovae, 
already constrain the deviation of w from -1 for flat mod- 
els better than to 10 % dSpergel et alj|2006h ISeliak et alj[2006t 
iTegmarke t al. 20061 and others.) Whether or not future obser- 
vations continue to converge on w = — 1, the dynamics of per- 
turbations will be crucial in elucidating the nature of cosmic 
acceleration. 

Even if w(z) = — 1 at low redshifts, this does not neces- 
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sarily imply a co smological constant. Certain models (e.g. 
Fard on et al.1 2006) predict that w = — 1 at the present epoch, yet 
the dynamics of perturbations differs at high redshifts. Contrary 
to common belief, it is even conceivable that w « — 1 yet low- 
redshift perturbations deviate from those in the ACDM model. 
On the other hand, if w 1 then perturbations of dark energy 
in the inhomogeneous metric are unavoidable. (Even if dark 
energy appears unperturbed in one spacetime slicing, its pertur- 
bations for w^-l are necessarily nonzero in any different slic- 
ing.) Exploring the perturbations' properties, specifically the 
properties considered in this paper, will then become pivotal to 
establishing the nature of dark energy. 

Density of dark radiation 

The current observations of CMB temp erature and polariza- 
tion, including th e WMAP 3-year results (Hins haw et al.| [2006; 
Page et al. 200 6)), when combined with the SPSS gala xy power 
spectrum (Tegmar k et al.ll2004t lEisenstein et alj|2005l) . with or 
without Ly-q fore st data, prefer enhanced neutrino density. 
ISeliak et alj d2006l) report N v = 5 .3^ at 2a, disfavoring the 
standard val ue N v = 3.04 at 2.4 a . The WMAP team in its lat- 
est analysis (Sper gel et al.l 1200 6) likewise concludes that from 
the WMAP3 and SDSS data N v = 7.1^ J. A similar prefer- 
ence for high N„ by the WMAP3 and S DSS data is seen by 
ICirelli & Strumial d2006l) . although not bv lHannestad & Ra ffelt 
(120061) . 

While neutrinos noticeably speed-up the background expan- 
sion in the radiation era, by its elf this leads to almost no observ- 
able cosmological signa tures dHu et al Jl^lBashinskT&S eli ak 
2004; Bashinsky 2004), given the fre edom of a compe nsating 
adjustment of matter density Q, m h 2 ( Bowenetal. 2002) and 



primordial helium fraction Y p (Bashinsky & SeliakT l2004l) 

This degeneracy is broken by the differences in the evolution 
of streaming neutrino perturbations and of the acoustic waves 
in the photon-baryon fluid. It is also broken by independent 
constraints on the helium fraction. Various signatures of neu- 
trinos that then appear remain partly degenerate with nuisance 
parameters, such as the dark energy equation of state or the 
normalization and tilt of the primordial power. Yet these signa- 
tures and the constraints on N v with the correspondingly differ- 
ent experiments have rathe r incomparable degree of robustness 
dBashinskv & Seliakll2004l) . 

The preference of WMAP3+SDSS data for an enhanced neu- 
trino density could be due to the physical excess of matter 
power over CMB p ower, expected for higher density of freely 
stream ing particles (Hu & Sugiva mall996l ;fBashinskv & Seljak 
2004). Yet other explanations, e.g., incomplete treatment of 
recombination or insufficient accuracy of bias modeling, also 
cannot be definitely excluded. 

In addition to increasing the ratio of matter to CMB power, 
freely stream ing neutrinos shift the p hase of CMB acoustic 
oscillations (Bashi nsky & Seliakll2004l) . The phase-shift sig- 
nature is very robust to systematic uncertainties and will be- 
come the primary discriminative mechanism for the near-future 
tight CMB constrai n ts on neutrino density dLopez et aTlfl 999 ; 

200^ 



Bowen et alj I2002t iBashinskv & Seliakl [2004; iPerotto et al.l 
2006), e.g. a(N v ) ~ 0.2-0.3 expected from the Planck mis- 
siorQ. It can provide decisive evidence for the excess of stream- 
ing relativistic particles or rule this possibility out. 

1 http://www.rssd.esa.int/Planck 



Nature of dark radiation 

If there is a true excess of energy density in the radiation 
era, at least three alternatives for non-standard dark radiation 
are possible: relic decoupled particles, self-interacting parti- 
cles dChacko et alj|2004 120051; beacom et alJl2004tlH annestad 



2005 [), or a tracking clas s ical field dRatra & Peebles! 119881; 
|p. Ferreira. M. JovcdTl997t IZlatev et al.l Il999t) . The pertur- 
bations of dark radiation propagate differently in all of these 
scenarios. As follows from this work, this in principle allows 
their experimental discrimination. 

1.2. Questions not answered by black-box computations 

For a particular cosmological model it is generally straight- 
forward to calculate linear power spectra and transfer functions 
with standard codes, if necessary, modified to include new dy- 
namics. Despite this, numerical calculations have limited use- 
fulness for exploring the signatures of new physics: First, for 
typical models, with close to ten unknown nuisance parame- 
ters, it is often intricate to establish numerically which of the 
observable signatures of the new physics cannot be compen- 
sated by parameter adjustments. While only such signatures are 
the true discriminators of the physics in question, they may be 
tiny and easily overlooked among large yet degenerate effects. 
Moreover, for every extended parameterization of the nuisance 
effects as well as for any additional constraining experiment, a 
new numerical analysis is required. 

Equally importantly, a numerical black-box computation 
does not answer which aspects of a model are responsible for 
its observable signatures. This obscures the separation of the 
physical facts that are backed by observations from the models' 
features that are believed true yet are untested. 

1.3. Approach, outline, and conventions 

The approach of this paper is to explicitly track the evolu- 
tion of gravitationally coupled inhomogeneities of the visible 
and dark species. This allows us to identify which observables 
are affected by which of the various properties of the dark ki- 
netics. Importantly, this approach also reveals the mechanisms 
behind the sensitivity of cosmological observables to various 
dark properties. Establishing the mechanism allows us to judge 
more intelligently the robustness of cosmological constraints, 
particularly, to know when this mechanism may fail or produce 
a different outcome. 

As we will see soon, there is an important subtlety in per- 
forming such identifications. To decipher the signatures of the 
dark kinetics, it is essential to address the gravitational interac- 
tion of perturbations during horizon entry. Then and only then 
perturbations of all abundant dark species are gravitationally 
imprinted on the visible species without suppression. The sup- 
pression of the impact of the species' energy overdensity Sp/p 
on the metric on subhorizon scales is evident from the Poisson 
equation 

$ = -k~ 2 A-KGa 2 8p ~ (H/k) 2 Sp/p, (1) 

where k is the (comoving) wavevector of a perturbation mode 
and H is the expansion rate (in conformal time). The impact of 
the species' velocities is suppressed even more [by (Ti/kf, c.f. 
eq. (fTTb.l At horizon entry, on the other hand, the factor H/k 
approaches unity and does not cause suppression. In particu- 
lar, only during the horizon entry are the visible species influ- 
enced by the perturbations of dark radiation and dark energy, 
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for which the Jeans length is expected to be close to the horizon 
size. 

Numerous authors (e.g., [ Press & Vishniac 1980; iBardeenl 
1980t iMa & Bertschingeril 19951) have pointed out considerable 
freedom in representing the inhomogeneous evolution on the 
scales of the order of and exceeding the Hubble scale. This 
freedom is due to, literally, an infinite number of possibilities 
for coordinate gauges, as well as for the variables (even the 
gauge-invariant ones) that can parameterize large-scale pertur- 
bations. In principle, this descriptional freedom could introduce 
ambiguity in relating observable features with specific physi- 
cal mechanisms that operate on large scales. Moreover, such 
relations indeed differ among the authors who use apparently 
dissimilar while formally equivalent descriptions — differ sub- 
stantially for some pronounced and important for constraints 
features. 

We argue in Sec. [2] that, as far as the observable impact of 
dark dynamics is concerned, there is little room for ambigui- 
ties. A well defined distinction can be drawn between an appar- 
ent connection "dynamical cause — > observable effect" that is a 
descriptional artifact or that is an objective causal relation. We 
will also see that in certain formalisms, which exist and can be 
distinguished by simple criteria, the physical microscopic prop- 
erties that characterize species at a particular time influence the 
apparent large-scale perturbations at the same time. Such for- 
malisms markedly simplify mapping of the characteristics of 
cosmological observables to the responsible microscopic dark 
kinetics. 

In the subsequent sections we do such a mapping for the 
CMB temperature power spectrum and LSS transfer functions. 
In Sec. [3]we discuss the general features of the primary probes 
of inhomogeneous cosmological dynamics. We also give dy- 
namical equations to be used further to quantify the CMB and 
matter response to dark parameters. In Sec. |4]we review various 
general parameterizations of potentially accessible information 
about the properties of dark sectors. We consider the parame- 
terizations of the metric, of dark densities and stresses, and of 
internal dark dynamics. Sec. [5] is central to our study. In this 
section we identify the characteristics of the CMB and cosmic 
structure that reveal such general properties of dark species as 
anisotropic stress, stiffness, clustering, and propagation of in- 
homogeneities. In Sec.|6]we study the specifics of cosmologies 
with modified gravity and discuss how modified gravity can be 
distinguished observationally. Our main results are summarized 
in Sec.|7]and its TableQ] 

Throughout the paper, distances will be measured in comov- 
ing units. Evolution will usually be described in conformal time 
dr = dt I a, where a is the cosmological scale factor and dt is 
proper background time. Overdots denote derivatives with re- 
spect to the conformal time, and Ti = a/a = Ha gives the Hub- 
ble expansion rate in conformal time. 

2. IS MAPPING OF LARGE-SCALE DYNAMICS TO 
OBSERVABLES UNAMBIGUOUS? 

Except for Sec. [6] we will consider models in which dark 
and visible sectors couple by standard Einstein gravity. These 
models correspond to a local action S = J d 4 Xy/=g£ with C = 
£dark+£ v i s +£ grav , where the dark degrees of freedom are repre- 
sented by the Lagrangian density £<jark, the visible ones by £ v ; s , 
and the gravitational ones by £ grav = (\6irG)~ l R. In agreement 
with observations, we assume that the large-scale evolution 
can be presented as a mildly perturbed Friedmann-Robertson- 
Walker (FRW) expansion. For this section only, we clock 



the evolution by the number of e-foldings N = \na, tractable 
through local matter density or CMB temperature. 

Let [N\,N2] be an evolution interval in the past. Consider a 
class, Mn, of all conceivable models with same £ v ; s that: 

(i) evolve identically for N <Nt', 

(ii) during the interval [N\ , N2] possibly differ in the micro- 
scopic laws of their dark, but not visible, dynamics; and 

(iii) by N = Ni have identical distribution of the dark species 
(in some non-degenerate measure) among all M\2 mod- 
els. 

By N = N%, the distributions of the visible species among these 
models will generally differ: The visible species would gen- 
erally be gravitationally affected earlier by the model-specific 
evolution in the dark sectors 

Given a quantity P that characterizes internal dark dynamics 
and an observable O, the following may be true: Regardless 
of the value of P during the interval ANn = [M,A^], all the 
models in the above class M\2 have the same observed value 
of O. Then it is natural to say that the observable O is not 
sensitive to the considered dynamical property P in the interval 
AN n - 

In linear order, this applies to the effects of any dark prop- 
erty P on any observable O that depends o nly on the perturba- 
tions that were superhorizon during ANn- IWands et alj (|2000) 
showed that perturbations in several gravitationally coupled flu- 
ids can be described by variables which during superhorizon 
evolution (for non-decayingmodes) are time-independent in 
each of the individual fluids|j This result can be extended be- 
yond multifluid models to any sector that is perturbed (inter- 
nally) adi abatically, i.e., is homog eneous in at least some coor- 
dinates dBashinskv & Se liak 2004). In other words, linear per- 
turbations encode no information about the microscopic prop- 
erties that characterized the dark universe when these perturba- 
tion were superhorizon. 

The abundances and kinetic properties that dark species have 
since horizon entry do influence the observed CMB or mat- 
ter perturbations. However, because of the aforementioned 
freedom in representing the evolution of large-scale inhomo- 
geneities, in many formalisms (including the most popular 
ones), the properties characterizing the dark species at horizon 
entry affect the apparent perturbations of the CMB or matter 
long before or after the entry. This can mislead one into viewing 
an observable feature as a probe of an entirely unrelated epoch 
and/or physical process. This source of existing and potential 
new misassignments can be eliminated systematically and nat- 
urally as follows. 

The impact of local dark dynamics on perturbations of visi- 
ble species will appear concurrent with the underlying micro- 
scopic dark ph ysics in any form alism which has the following 
two properties (Bashinsky 20061): 

2 Even the dark species in the M\i models may not be perturbed identically 
at Ni in terms of every measure. On large scales the quantification of dark 
perturbations by a different measure can depend on the perturbations of visible 
species and the metric. We require only that the dark perturbations are the same 
in at least one non-degenerate measure. 

3 Linde & Mukhanov (2006) noted that under certain conditions, which may 
exist, e.g., during reheating, gravitational decays of species into another type 
of species may mix the superhorizon perturbations in different sectors. Such 
decays are, nevertheless, negligible for the post-BBN evolution studied in this 
paper. 
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I. Perturbations are frozen on superhorizon scales. 

II. Perturbations evolve by the equations of the FRW metric 
whenever the geometry is unperturbed. 

We imply that the description of gravity in these formalisms re- 
duces to the Newtonian one in the weakly perturbed metric on 
subhorizon scales. This is easily achieved, e.g., by parameteriz- 
ing metric inhomog eneities by the gravitational potentials of the 
Newt onian gauge (Mu khanov et al.|[l992t iMa & Bertschi nger 
fl995h . 

Then the apparent linear impact of dark species on observ- 
ables will be found to be identical in any description with prop- 
erties I and II: This impact is practically unambiguous for the 
evolution of perturbation modes after horizon entry; the appar- 
ent perturbations do not evolve at all prior to the entry by con- 
dition I. Therefore, these descriptions could differ only by the 
changes of pert urbations during horizon entry. It is straightfor- 
ward to argue ( Bashinsky| |2006l) that there are no such differ- 
ences in linear theory among the formalisms that satisfy I and II. 

One such particular, natural and simple, formulation of the 
full linearly perturbed Einst ein-Bol t zmann cosmological dy- 
namics was described in Bashinsky (2006). This formulation 
is based on canonical coordinate variables. The canonical vari- 
ables for perturbations of phase-space dis tributions or radiation 
intensity wer e considered in the past (e.g., Chandrasekhar 1960; 
iDurreri 120011) but have not become mainstream in the present 
cosmology. In addition to the demonstrated more direct con- 
nection of large-scale evolution to microscopic kinetics of cos- 
mological species, this formulation has several technical advan- 
tages over popular alternatives, most of which consider pertur- 
bations of proper quantities. It is used in the analysis that fol- 
lows. 

3. PROBES 

Following the standard route, we separate the effects of 
dark sectors on the metric into the contributions to the back- 
ground geometry and to metric inhomogeneities. The back- 
ground geometry is fully described by the Hubble expansion 
rate as a function of redshift and by possible global curvature. 
It is constrained by a variety of "geometrical" probes, e.g., the 
luminosity-redshift relation for type la supernovae or the an- 
gular size of the acoustic features in the power spectra of the 
CMB and matter. In this paper we focus on extracting the po- 
tentially much richer information about the dark dynamics that 
is contained in the metric inhomogeneities. This information is 
to be inferred from the imprints of the metric inhomogeneities 
on such various observables as the CMB, galaxy distributions, 
lensing shear, quasar absorption spectra, etc. 

The primary probes of dark perturbations can be broadly 
classified as either "light" or "matter": The "light" probing 
the metric along trajectories which are close to null geodesies 
(ds 2 ~ 0, e.g., CMB photons); the "matter" moving almost 
with the Hubble flow with non-relativistic peculiar velocities 
(\dx/dr\ -C 1, e.g., galaxies orLy-a clouds). These two classes 
of probes, considered in Secs. |3.1| and |3.2| respectivelv, provide 
information that is complimentary in several respects, discussed 
in details in subsequent sections. 

Linear theory adequately describes the inhomogeneous dy- 
namics during horizon entry, when physical quantities are 
perturbed by about 0.001%, and long after the entry. Then 
the signatures of the dark kinetics are encoded in the lin- 
ear transfer functions [alternatively, in Green's functions 



(Ma gueiiol 119921: fe accigalupi 1999; iBaccigalupi & Perrottal 
20001: iBashinskv & Bertschingerl 1200 ll |2002|) 1, which are the 
essential constituents of "dynamical" cosmological observ- 
ables, such as power spectra, luminosity functions, etc. To 
be specific (though general within the linear regime), in this 
section we consider an individual perturbation mode with a 
comoving wavevector k. 

3.1. CMB 

After horizon entry until photon decoupling around hydro- 
gen recombination, a perturbation mode in the photon-baryon 
plasma undergoes acoustic oscillations. All the memory of the 
gravitational impact at the mode's entry is then retained only in 
its oscillation amplitude A(k) and its temporal phase shift (f(k). 
These quantities map into the heights and phase of the acoustic 
peaks in the observable CMB angular power spectra Q, with a 
rough l^rk correspondence I ~ kr^, where rd ec is the angu- 
lar diameter distance to the CMB decoupling surface. [In the 
standard ACDM model r dec = r(z dec ~4 100) ~ 14Gpc] 

Our goal is to establish how A(k) and ip(k) are affected by 
the gravitational impact of perturbations in various species. As 
noted in the introduction, this task involves a subtlety that with- 
out proper care can lead to erroneous conclusions. Although 
the notion of density or temperature perturbation is uniquely 
defined in the FRW geometry, ambiguities arise on large scales 
in the perturbed metric. We argued in Sec. |2]that this compli- 
cation is uniquely resolved in linear theory, where switching 
off a microscopic effect of interest results in the same change 
of observables regardless of other local properties of the vari- 
ous species. We can then calculate the observables in the mod- 
els with the effect "on" and "off" and identify the signature of 
the effect with the difference. Being interested in the kinetics 
of perturbations, we compare the models with identical back- 
ground expansion H{z). For the compared models, we assume 
identical initial conditions, i.e., identical dynamics during the 
inflationary epoch, when the superhorizon perturbations likely 
were generated. 

We d escribe the perturb ative evolution by a formalism de- 
tailed in Bashinsky] (120061) . Then the gravitational forcing of 
perturbations appears concurrent with the responsible local in- 
teractions and the above unambiguous causal relations are man- 
ifest. In this formalism the general-relativistic generalization of 
particle overdensity is d = Sn coo /n coo , a perturbation of particle 
number density per coordinate volume. 

When the CMB is tightly coupled to electrons by Thompson 
scattering, its (coordinate) overdensity evolves as 

d\+[^- b H + 2T d k 2 ^d 1 + c 2 k 2 {d 1 -D) = Q (2) 

(Bas hinsky & Seljald[2004h . Here and below R b = 3p b /(4p 1 ), 
the diffusion t ime T £ / = [l-||(l+^ i )" 1 +(l + Rh T 2 ] /(6aw,.O Thomp) 
dKaisedll983"l) determines the Silk damping dSilkH T968). and 
c 2 = [3(1+ Rb)]~ l gives the sound speed in the photon-baryon 
plasma. Only scalar perturbations are considered^ The 
gravitational driving of the CMB modes in eq. (f2]i is medi- 
ated through an instantaneous equilibr ium value D of photon 
overdensity. In the Newtonian gauge (Mukha nov et al J fl992h 

4 The connection between dark dynamics and its observable signatures is rela- 
tively straightforward in the tensor secto r, where gauge ambiguities are absent 
i Bardeen 1980; Kodama & Sasaki 1984). Vector perturbations, even if primor- 
dially generated, are not expec ted to survive superhorizon evolution (Bardeen 
1980; Kodama & Sasaki 1984). If necessary, they can be analyzed similarly to 
the scalar perturbations. 
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iMa & Bertschingeril 19951) 

ds 2 = a 2 [-( 1 + 2<5>)dT 2 + ( 1 - 2^)dx 2 } , (3) 

the driving term equals 

D(T,k) = -3($ + V+R b ®). (4) 

In the radiation era, when SjCl, the driving of the CMB by 
dark perturbations on all scales is controlled only by the sum 
This also applies to the CMB photons on all scales after 
their decoupling from baryo ns, c.f. the equ ation for the evolu- 
tion of CMB intensity ©. iDurreri d 19941) argued that this is 
natural for the relativistic particles, whose dynamics is confor- 
mally invariant, and therefore should be sensitive only to the 
Weyl part of the curvature tensor. Inde ed, the scalar part the 
Weyl tensor is fully specified by $ + \l> dDurrerill994l) . In the 
decoupled case, the driving by $ + ^> is well known as the in- 
tegrat ed Sachs-Wolfe effect (ISW) (|R. K. Sachs. A. M. Wolfd 
fl967h . 

The full linear d ynamics of partially polar ized CMB photons 
and baryons (e.g., Ma & Bertschingerf l995) is also straightfor- 
ward to formulate so that the changes of the apparent inhomo- 
geneit ies are concurren t with the responsible local dark prop- 
erties (Bashinsky 2006). In particular, the perturbation of the 
polarization-summed CMB intensity corresponds to its canon- 
ical variable Ljx^. rij) = /~/Zy — 1 that evolves a ccording to a 
transport equation ( Dur rerl 1 990t lBashinsky||2005l also present- 
ing a fully nonlinear treatment) 

i+ mVn = -4«,V,($ + *) + C T , (5) 

where Cj is the Thompson collision term. 

3.2. Matter 

The evolution of linear scalar perturbations of cold dark mat- 
ter (CDM) on all scales and, after decoupling from the CMB, 
of baryons on large scales is governed by the conservation and 
Euler equations 

d c + d iV ' c = 0, vi+Hv^-di®. (6) 

Their solution for a Fourier mode with non-singular (inflation- 
ary) initial conditions is 

2 r 1 r 

d c = d cm — k / u c dr, u c =- (a^)dr. (7) 
Jo a J 

The quantity u c that appears in these equations is physically the 
CDM velocity potential: v' c = -diu c . 

Similarly to the CMB modes, the matter perturbations carry 
the memory of the metric inhomogeneities at the horizon entry 
in two independent functions: the overdensity change and ve- 
locity boost that were generated during the entry. After horizon 
entry, the excess of CDM velocity, decaying as 1 /a, yet not van- 
ishing instantaneously, continues to enhance matter clumping. 
Hence, the implications of the velocity boost for the observed 
5pm/ Pm are genera lly more important tha n the initial overden- 
sity enhancement dHu & Sugivamafll996l) . 

Eqs. (0 show that matter p erturbations are driven by t he po- 
tential $ alone on all scales dBashinskv & Seriakll2004l) . This 
also has a natural explanation. In the Newtonian gauge, the 
point particles of a mass m are described by the Lagrangian L = 
-mV-ds 2 /dr = -ma^[l+2$(x)]-[l-24'(x)]x 2 . The corre- 
sponding contribution of '5 is negligible for nonrelativistic par- 
ticles, when L ~ ma [x 2 /2-<I>(x)] + const. Moreover, the initial 
(superhorizon) perturbations of particle's coordinates x deter- 
mine the initial value of the coordinate overdensity d c also in- 
dependently of 'J. 



Since the CMB, on the other hand, is affected by the 
sum <f> + "J, this is one of the respects in which matter and the 
CMB provide complimentary information about dark inhomo- 
geneities. A combined analysis of both probes is particularly 
useful for testing a deviation of \] //<i> from unity, produced 



by freely streaming species (e.g., IMa & B ertschinger 



1995) 



2006; 



or generally expected in m odified gravity dBertschingeri 
iDodelson & Liguorill2006l) . 

4. PARAMETERIZING DARK DYNAMICS 

In this section we review several general parameterizations 
of the potentially observable properties of the dark universe. 
We start from the spacetime metric, connected most directly to 
observational probes. Then, assuming validity of the Einstein 
equations, we match the metric characteristics to a general de- 
scription of dark densities and stresses. Finally, we consider 
a parameterization of the potentially measurable properties of 
dark internal dynamics. The observational probes of these pa- 
rameters will be studied in Sec. [5] The status of these param- 
eterizations in cosmologies with modified gravity will be ad- 
dressed in Sec. [6] 

4.1. Metric 

Provided the visible matter couples covariantly and mini- 
mally to a certain metric tensor g M „, the observable impact 
of dark sectors or modified gravity alike can be parameterized 
by a single function of redshift z, to describe the FRW back- 
ground, and by two independent functions of z and k, to specify 
scalar metric perturbations. For example, the background can 
be quantified by its expansion rate H(z) (with its possible spa- 
tial curvature), and the metric perturbations by the Newtonian- 
gauge potentials $(z,k) and *(z,k). 

4.2. Densities and stresses 

For a more direct link to the internal properties of the dark 
sectors, it may be worthwhile to parameterize their contribu- 
tion to the energy-momentum tensor = J2 a ^a V ■ ^ gravity 
is described by general relativity then determines the geo- 
metrical quantities {H, $, \I r } through the Einstein equations. 
In this view, the background may be specified by the average 
energy densities of the species p a (z) = ~^o(ay ^he scaiar inho- 
mogeneities of the metric are determined from the perturbations 
of the energy densities, sourcing the curvature potential "J, and 
from anisotropic stress, generating the splitting $-5'. The cor- 
respond ing formul as are given next. 

After |Hu|C998), we will use proper particle number over- 
density (5^ c) (z,k) = Sp^ /(p„ + p a ) in the comoving gaug^fl In 
any other gauge 



5 f= JP^+3Hu a , 

Pa+Pa 



(8) 



where v' a = -diu a . Specifically, 8^ can be expressed in terms of 
the dynamical coordinate overdensities d a , used in Sec. [3] as 

S^=d a + 3Hu a +3^. (9) 

The comoving-gauge overdensity 5^ is convenient for two 
reasons. First, it directly sources the Newtonian potential ^ 
dBardeenll 19801) 

V 2 ^ = ^Ga 2 Y J (p a +Pa)8^\ (10) 

a 

5 The comoving gauge conditions for scalar perturbations are zero total mo- 
mentum density (7j° = 0) and vanishing shear of the spatial metric (gij oc 
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where V 2 = (3) g' ; <9,<9; is the spatial Laplacian (from now on the 
background is assumed flat). Second, the generalized pressure 
gradient in the Euler equation [the term -d-J p in eq. ( fL9j >l in the 
comoving gauge takes its usual special-relativistic form [f p = 
Sp^/ip+p)]. 

Nevertheless, we view <5jf' only as a useful linear combina- 
tion (O, rather than an independent dynamical variable: The 
evolution of S ( f\ unlike that of d a , is highly nonlocal. Beyond 
the horizon S ( f\ unlike the frozen d a , is affected by all species, 
in addition to the studied dark species a. In terms of the dynam- 
ical variables d a and u a the Poisson equation ( fTQb reads 



V 2 \ V 

3 + — Ur = 2_^x a (da + 3Hu a ), 



(11) 



with 



7 = 4nGa\p+p) = - (1 + w)W 



, Pa + Pa 

and x a = . (12) 

p+p 



Unlike eq. ( fTQb . the generalized Poisson equation in the 
form dTTb in the superhorizon limit is explicitly non-singular 
and its solution W -^(d+3TLu) even becomes local. 

We describe the scalar component of the anisotropic stress 



Tj - by a potential a(z, k) as 



(13) 



3 (Pa+Pa) 

The corresponding gravitational equation is dBardeenl [l980; 



iKodama & Sasakifl984tlMa & Bertschingerill995l) : 
$ = * 



(14) 



Thus the densities and stresses of the background and scalar 
perturbations may be parameterized by a set {p a (z), (5jf'(z,k), 
tr a (z,k)} for each sector a. Both 5^ and a a , as defined by 
eqs. (O and ( fT3l ). are gauge-invariant quantities. Neither of 
them, though, can be measured locally. 

4.3. Dynamical properties 

As a step further, we can try to parameterize the dynamics of 
the dark densities and stresses. One such u seful a nd popular pa- 
rameterization has been suggested bv lHul (119981) . He described 
the dynamics of the background by the conventional "equation 
of state" w a (z) = p a /Pa, and the dynamics of perturbations by 
the effective sound speed 

c^ a {z,k) = 8pf{k)/8pf{k). (15) 

We call this quantity "stiffness" throughout the paper in order 
to distinguish it from the physical velocity of perturbation prop- 
agation c p , Sec. 15.31 which may differ from c e s and is probed 
by different observables. 

In addition, iHul (1998) considered a "viscosity parameter" 
c 2 is a (z,k) (defined by the ratio of the dynamical change of 
anisotropic stress to the velocity gradient). However, realisti- 
cally the evolution of anisotropic stress is determined not only 
by the velocity gradient but also by additional degrees of free- 
dom, e.g., by the I = 3 multipole or polarization. Therefore, 
in our study we will not use c 2 is but will work directly with 
a a (z,k). 

A dark sector that couples negligibly to the other species can 
be assigned a covariantly conserved energy-momentum tensor 
Tq_ through the s tandard procedure of metic variation (e.g., 
Land au & Lifshitzll 19751) . Given the value of the sector's back- 
ground energy density p a at some moment (e.g. the present 



density, set by VL a h 2 ), w a (z) determines the density at all times 
from 

p = -3H(l + w)p (16) 
(the subscript a is dropped from now on.) The evolution of den- 
sity perturbation d follows by linearly perturbing the equation 
of energy conservation T 0l/ . v = 0: 

d + dy = Q (17) 

witrQ 

p6p-pSp_ 



(p+pY 



-3H 



eff 
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(18) 



Linearization of T w .„ = yields the evolution of momentum- 
averaged velocity^ 

v i +nv i =-d i ^-d i f p +F i a (19) 

with 

Sp-pu 

Jp 



P + P 



- eff 



(20) 



P+P 3 



(21) 



We remember that <5 (c) is given by eq. (0, and the gravita- 
tional potentials are determined by the linearized Einstein equa- 
tions ( [Tol l and (TPfl) . These equations are closed after c 2 ff and a 
are specified by the species' internal dynamics. The expansion 
of perturbations over Fourier modes splits the equations into 
uncoupled systems for each mode, with its own c 2 ff (A:) and a(k). 

Similarly to cr{k), the stiffness c 2 ff (£) is manifestly gauge- 
invariant but not measurable locally. 

4.4. Other parameterizations 

Any of the above sets {H, {p, S (c) , a}, or {w, c 2 ff , a} 

parameterizes all properties of arbitrary dark sectors that can 
be constrained by any probe of background and scalar pertur- 
bations. It may nevertheless be useful to explore other internal 
properties of the dark species. Certain properties, e.g. particle 
masses and cross-sections, that are not readily extractable from 
the above sets may be important for particle-physics models. 
Other properties may have more distinctive or less degenerate 
observable signatures. An example appears in Sec. 15.31 dis- 
cussing the dark characteristics that control the additive shift of 
the CMB peaks. 

Of course, yet alternative parameterizations of the i nternal 
dark dynamics may also be considered. For example, iLinderl 
(2005) suggested to approximately quantify the growth of cos- 
mic structure by a growth index parameter 7. This approach 
may be useful for utilizing LSS to probe modified gravity or 
non-standard long-range interaction. We do not pursue it here, 
being interested in more universal descriptions that are applica- 
ble to the CMB as well. 

5. MAPPING DARK DYNAMICS TO OBSERVABLE FEATURES 

In this, main, section we identify the observable signatures 
of the dark properties. Specifically, we consider the signatures 
of: anisotropic stress, the dark species' stiffness, the propaga- 
tion speed of their inhomogeneities, and their clustering. The 
observational probes considered are the transfer functions for 
cosmic structure and the power spectra of the CMB. 

p/p, appearing in eq. I18K can be related to the background equation of state 
was p/p^w-^^. 

1 Eqs. 4 1 9i — 121 ) are simpler in terms of the velocity potential: for v' = — d\u, 
they read 

ii + Hu = <S> + c£ ff 8 {c) - - V 2 ct. 
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FIG. 1 . — The evolution of CMB overdensity dj(r, k) (oscillating, red curves) and the corresponding gravitational driving term D = -3(<I > + *&) (monotonically 
falling, brown curves) in the radiation era for Rf, <C 1 and negligible photon diffusion. The gravitational driving of the CMB [eq. (2)] is illustrated with a spring that 
connects the d 7 and D curves. The panels demonstrate the impact of perturbations in decoupled neutrinos (left) and quintessence (right). — Left panel shows d 1 
and D for: tightly-coupled photons only (dashed), the standard model with three neutrinos (wider solid), and the same model after switching off neutrino anisotropic 
stress as a source of <3?-'J' (thinner solid). — Right panel displays: the earlier photon-only model (dashed) and a fictitious model where, as in the standard model, 
59% of energy density is in photons but the remaining 41% is now in a tracking quintessence </>, with wj, = 1/3 (solid). 



5.1. Anisotropic stress 

Already in the standard cosmological model, significant 
anisotropic stress is generated in the radiation era by freely 
streaming neutrinos. The anisotropic stress in our uni- 
verse may differ from the standard three-neutrino expec- 
tation because of, for example, additional free-streaming 
relativistic species, which would enhance the stress. Or, 
it may differ because of non-minimal neutrino cou plings 
(IChacko et all |2004 l2005t iBeacom et alJ 12004 lOkuil [2005; 
iHannestadl 120051 : iGrossman et al.ll2005l) . which would locally 
isotropize neutrino velocities and reduce or eliminate their 
anisotropic stress. Aniso tropic stress a changes $ and 'f 
dMa & Bertschingerlll995l) and consequently the CMB gravi- 
tational driving term D (01 even on superhorizon scales. Thus 
the a of dark species affects the CMB modes since the earliest 
stages of horizon entry. 

The left panel of Fig.Q]shows the evolution of gravitationally 
coupled photon and neutrino perturbations in the radiation era 
for various assumptions about neutrino abundance and prop- 
erties. The plots are obtaine d by numer i cally integrating the 
corresponding equations from lBashins kv (2006). The oscillat- 
ing, red curves give the CMB overdensity d 7 (T,k), The falling, 
brown curves show the driving term D(r,k), equal to — 3(<E» -I- M/ ) 
in the radiation era. 

As evident from this figure, the photon overdensity d~ in the 
presence of freely streaming neutrinos (wider solid oscillating 
curve) is suppressed with respect to d 1 of a neutrinoless model 
(dashed oscillating curve). The suppression of the CMB pertur- 
bations by neutrinos can be related to the reduction of D early 
during horizon entry, with the CMB modes then experiencing a 
smaller initial boost. 

The suppression of D and, consequently, of the CMB os- 
cillations in the Newtonian gauge is due to er„ being a direct 
source of gravitational potentials. It is not due to the damping 
of neutrino perturbations by the viscosity effect of a v , eqs. ( TT~9b 
and d2lT i. We can verify this by repeating the calculations for 
the standard model with three neutrinos but with the following 
modification: We remove the gravitational effect of anisotropic 
stress in eq. dl4l >. i.e., set $ = \&. We make no other changes in 



the integrated equations, i.e., preserve the Poisson equation ( flOb 
and use the standard dynamical equations for all species' per- 
turbations. In particular, we retain the a v viscosity term in the 
v v equation [c.f. eqs. ( fl9ll2lT il. The results of this calculation, 
with the "gravity of a v removed'Q, are plotted on the left panel 
of Fig. Q] with the thinner solid curves. Switching off the direct 
gravitational effect of er„ is seen to lift the neutrino suppression 
of the subhorizon CMB oscillations. 

On the scales that enter the horizon in the radiation era 
(£ ^> 200) the suppression of the amplitude of the CMB os- 
cillations, Ary, can be cal culated analytically in le ading or- 
der in R v = Pu/(p~ f + Pv) dBashinskv & Seliakll2004) . For all 
physical values < R„ < 1, this suppression is well fitted by 
( lHu&Sugivamalll996i) 

s^r( 1+ Kf' (22) 

where R v is vari ed at unchanged flu ctuations of the primor- 
dial curvature £ dBardeen et al.l fT983). This variation implies 
fixed inflationary physics, although then the superhorizon val- 
ues of both $ and ^ change with R v . The oscillations in the 
CMB power spectra of temperature, polarization, or their cross- 
correlation for £ ^> 200 are suppressed by the square of d22b . 

To summarize, the suppression of the CMB oscillations by 
freely streaming neutrinos is due to neutrino anisotropic stress. 
In the Newtonian gauge the suppression may be attributed to 
the anisotropic stress acting as a source of metric perturbations, 
impacting CMB photons. The damping of bulk velocity and 
density perturbations of neutrinos by streaming plays almost 
no role in the suppression of CMB fluctuations. 



8 The performed removal of cr„ from the gravitational equation )14t but not 
from the Euler equation U9t violates the covariance of the Einstein equations. 
Therefore, the separation of the effects of anisotropic stress into "gravitational" 
and "dynamical" should not be considered beyond the context of the Newtonian 
gaug e. W e co uld preserve the Einstein equations by removing cr„ from both 
eqs. {14} and U9t . The result would be the neutrinoless model, described on 
the left panel of Fig.[T]by the dashed curves; again, showing no suppression of 
the CMB amplitude. 
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5.2. Stiffness c 2 s 

The quantity c^ s (z,k) = Sp^'/Sp^ describes the stiffness of 
the dark medium to perturbations with a wavevector k. (We 
avoid calling c^ s "the sound speed," not to mix it with c p of 
Sec. 15.31 ) The stiffness c\ {{ affects directly the evolution of the 
dark species' overdensities and peculiar velocities. The direct 
impact of c 2 n on the evolution of overdensity, eq. <fl7l l. is de- 
scribed by the term Q ( fT8l . proportional to the so called "non- 
adiabatic" pressure Sp-(p/ p)Sp. The velocities are accelerated 
by pressure gradient, the term —dif p = -di(c 2 n 5 p^) in eq. ( fT9b . 
directly proportional to c 2 ^. For perfect fluids and k^S> Ti, the 
quantity c e g{k) gives the (phase) velocity of acoustic waves. In 
general, however, c 2 s need not be related to the speed of pertur- 
bation propagation, which will be studied in Sec. 15. 31 

The dependence of the dark perturbations' evolution on their 
stiffness c 2 n is reflected in the perturbations' contribution to 
metric inhomogeneities. Throug h the latter, the dark stiffness 
affects the visible species (e.g.. lErickson et a"Dl2002h . To be 
specific, we assume that superhorizon perturbations are adia- 
batic0 Then the comoving overdensity S (c \ eq. ©, vanishes 
in the superhorizon limit. Consequently, d a [eqs. ( 11711181 )1 and 
v'fl [eqs. ( TT9ll2"0l >l depend on c 2 s only at the order 0(k 2 T 2 ). The 
sourced potentials W and $ therefore become sensitive to c 2 n 
only at a late period of the horizon entry, also only at the or- 
der 0(k 2 T 2 ). Note that, as discussed in Sec. 15.11 anisotropic 
stress, in contrast, changes the potentials even in the superhori- 
zon limit kr — ► 0. Thus c 2 n affects the observable species at a 
noticeably later stage of the horizon entry than a does. 

The lateness of the impact of c 2 n is clearly seen on the right 
panel of Fig. Q] which shows the joint evolution of perturbations 
in a photon fluid and a classical scalar field <p (quintessence). 
The quintessence background density is set to track radiation 
(w0 =1/3) but its stiffness exceeds that of the photon fluid 

( c ls,<f, = 1 > 4r, 7 = !/ 3 > e.g.,E3(ll998|)). Either the pertur- 
bations of the quintessence (right panel) or of freely stream- 
ing neutrinos (left panel and previous Sec. 15.11 ) decrease the 
CMB-driving sum $ + \I>. With neutrinos, $ + '5 was decreased 
already on superhorizon scales by neutrino anisotropic stress 
(er„ 7^0 = cr 7 i). On the other hand, in agreement with the previ- 
ous arguments, $ + '5 is seen to be affected by the stiff inhomo- 
geneities of quintessence, whose anisotropic stress is zero, at a 
later time. 

There are several observable consequences of the lateness of 
the influence of dark stiffness on potentials. First, while the 
early reduction of <£> + 'J? by cr„ suppressed the amplitude of the 
acoustic oscillations in d 1 , the later reduction of $ + $ by large 
c eff <t> slightly enhances this amplitude. The right panel of Fig. Q] 
explains this paradox: The reduction of $ + ^ by quintessence 
perturbations increases the negative driving of photon overden- 
sity when d 1 is on average negative. As a result, the acoustic 
oscillations gain energy and their amplitude increases. 

The lateness of the stiffness impact has another important 
consequence. Namely, the ratio of the matter response over 
the CMB response is considerably larger to probing Cg ff than 
to probing a. This is evident from comparing the left panel of 
Fig. [2 showing CDM overdensity (black, rising curves), and 
the earlier Fig. [T] for the CMB overdensity, all in the radi- 

9 The following arguments remain valid under a considerably milder condition: 
It is sufficient that t he primordial perturbations in the probed dark species are 
internally adiabatic iBashinsky & Seljak 2004). In particular, this condition is 
satisfied automatically for any one-component fluid. 



ation era. In Fig. |2j CDM overdensity d c is plotted for: the 
photon-only model, where a = (dashed); the standard model 
with 59% of energy density in coupled photons and 41% in 
streaming neutrinos (wide solid); and a fictitious model with the 
same 59% of energy density in photons but the remaining 41% 
carried by quintessence (thin solid). We see that the enhanced 
stiffness of the dark component of the third model impacts the 
matter an order of magnitude more strongly than the anisotropic 
stress of the second model. On the other hand, the CMB per- 
turbations, studied on the left and right panels of Fig. Q] are 
affected by neutrinos and quintessence of the second and third 
models comparably. 

These results have a simple explanation. The early impact 
of a on matter velocities is washed out over time by the Hub- 
ble friction [Tiv' c term in eq. ©]. The Hubble friction, acting 
on massive CDM and baryons, damps less the c^ ff -dependent 
late impact of dark perturbations on matter. The overdensity of 
relativistic photons, on the other hand, is unchanged by Hubble 
friction [see eq. @, where Rt, is negligible]. Hence the present 
CMB anisotropy responds comparably to either of the impacts. 

The relatively high ratio of the matter over the CMB re- 
sponses to the dark stiffness should help distinguish an excess 
of relic relativistic particles in the radiation era from a sub- 
dominant tracking quintessence. Both the decoupled relativis- 
tic relics and quintessence have a similar n ondegenerate phase- 
shift signature in the CMB power spectra (Bashinsky & Seliak 
2004, more in Sec. 5.3 below). This signature will soon 
considerably tighten the constraints on neutrinos or an early 
quintessence with the enhanced angular resolution of CMB ex- 
periments. Yet it does not discriminate between the two sce- 
narios if a nonstandard signal is observed. The discrimination 
can, however, be achieved by combining the CMB data with ac- 
curate measurements of matter power, together sensitive to the 
strong effect of quintessence on the ratio of the CMB and mat- 
ter power. The discrimination should be facilitated by the fact 
that tracking quintessence and an excess of relativistic species 
change this ratio in opposite directions. 

5.3. The speed of sound or streaming 

Dark inhomogeneities were seen to affect the amplitude of 
the oscillations in the CMB temperature and polarization power 
spectra. In addition to the heights of the acoustic peaks, the 
CMB spectra are characterized by the peaks' positions £„. The 
period, of the acoustic oscillations is well known to be 
controlled entirely by background geometry and the proper- 
ties of the photon-baryon plasma. (Namely, £a = irrfe c /Sfe c , 
where ra ec is the distance to the surface of CMB decoupling 
and Sdec = Jq^ Cyb dr is the corresponding size of the acous- 
tic horizon.) On the other hand, the overall additive shift of 
the peak's sequence, i.e., the parameter A£ in an approximate 
relation £„ ~ uIa + A£ is a robust characteristic of dark pertur- 
bationsOB 

The shift of the peaks in £ is determined by the shift 
Aip of the temporal phase of the subhorizon modes, d 7 = 
A 1 cos(kc 1 t,T - Aip), as A£ = £aA(p/tt. This shift of phase is 
tight ly connected to the locali t y of the inhomog eneous dynam- 
ics dBashinskv & SeliakH2004tlBashinsk"vll2004l) . This becomes 
manife st in the real-space description of pertu rbation evolu- 
tion of IBashinsky & Bertschingeij d200U |2002|) . using plane- 



10 The positions of individual peaks receive additional corrections from the 
changes of the peaks' shape, such as due to the Silk damping and (for temper- 
ature) falling Doppler contribution. These corrections are fixed by the back- 
ground properties and should not spoil the constraints on the overall shift A£. 
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parallel Green's functions. The Green's function of J 7 in the 
radiation era, when the mode s forming the acoustic peak s en- 
ter the horizon, was found in iBashinskv & Sehakl (12001 . Its 
Fourier transform yields 

-W3 ($ + *)ui= c , T 



Aip'. 



(23) 



(for \A(p\ Here, the initial conditions are assumed adia- 

batic, and c, = c 7 b = 1 / V3 is the speed of sound in the radiation 
era. The term (^ + ^>)i x i =CsT is the value of <E> + <J' at the acous- 
tic horizon of a perturbation that is initially localized in space 
to a thin plane: d(r — > 0,x) = Sd(x), where 5d(x) is the Dirac 
delta function. For adiabatic initial conditions, ( < & + v I / )| A |=c,r 
is fully determined by the perturbations that propagate beyond 
the acoustic horizon. In particular, (^ + ^)\ x \=c s t = if none of 
the dark species support perturbations that propagate faster than 
the acoustic speed c s . [For a proof, accounting for subtleties of 
inhomogen eous gravitational dynamics on large scales, see Ap- 
pendix B of|BasSnskxASeli^(|200l)]. 

Two commonly considered types of cosmological species do 
support perturbations that propagate faster than c s . These are 
decoupled relativistic neutrinos and quintessence, in both of 
which perturbations propagate at the speed of light. The in- 
crease of the energy density in either of these species by one 
effective fermionic degree of freedom displaces the peaks of the 
CMB temperat ure and E-polarization spec tra by A/+i„ ~ -3.4 
for neutrinos dBashinskv & Seliakl 120 04) and Al + i v ~ -11 
for tracking quintessence dBashinskv! 120061) . This nondegen- 
erate effect will le a d to tight constrain t s on the abundance 
dLopezetalJ 1 1 999b lBowen etal.1 120021; feashinskv & Seliak 
2004 iPerotto et alj|2006l) and interaction Chacko et al.ll2004l) 
of the relic neutrinos from upcoming CMB experiments. These 
constraints notably improve with higher angular resolution in 
temperature and polarization as new data at higher ts allow 
better identification of the oscillations' phase. [The polariza- 
tion channel, whe re the acoustic peaks are m ore pronounced, is 
especially useful dBashinskv & Seliakll2004l) .l 

Thus a crucial characteristic that affects the acoustic phase is 
the velocity c p of the wavefront of a localized perturbation of 
dark species, shifting the phase if and only if c p > c s » 1 /v3- 
The velocity c p of a perfect fluid or a classical field is deter- 
mined by the stiffness c\ e . However, in general c p and c e ff are 
independent. As a notable example, free-streaming relativistic 
particles have c e ff(k) = 1 / V3 but c p = 1 . Unlike c e s(k), c p is de- 
fined irrespectively of any gauge and is locally measurable. Fi- 
nally, as discussed in Sec. l5.2l and the present section, the quan- 
tities c e ff and c p map to different observable signatures. Thus 
for gaining robust knowledge of the nature and kinetics of dark 
radiation, including neutrinos and possible early quintessence, 
it is important to constrain observationally both c e s(k) and c p . 

5.4. $ and ty; clustering in the dark sectors 

In the three previous subsections we analyzed the gravity- 
mediated signatures of the dark stresses. Now we shift our at- 
tention to the observational manifestations of potentials $ and 
'J themselves, irrespective of their cause or even of the valid- 
ity of the Einstein equations. As described quantitatively in 
this subsection, the impact of the gravitational potentials on 
large-scale structure (LSS) and the CMB differs in several es- 
sential ways. Consequently, the corresponding constraints will 
be rather complimentary. 

Before considering the observable impact of the metric po- 
tentials, let us summarize the scenarios in which the potentials 



are expected to differ from the predictions of the standard cos- 
mological model. As compared to a typical model with CDM 
and quintessence dark energy (including ACDM model as a 
limiting case), different $ and $ in a given background w(z) 
are expected: 

(a) During horizon entry for practically any other dark en- 
ergy model. 

(b) On any scales for models with modified gravity. 

(c) 



On subhorizon scales fo r models with non-canonical ki - 
netic term for the field ( Armendar iz-Picon et al. 200 lb . 
models with warm da rk matter dBlumenthal et al.ll 1 982t 

lOlive & Turnei 1982b or with interacting dark matter 

dCarlson et aljl992tlde Laix et alJl995l:ISpereel & Steinhardtl 
2000). More mundane physics affecting the subhorizon 
potentials includes thermal, radiative, or magnetic pres- 
sure on baryons, and various astrophysical feedback 
mechanisms (supernovae, central black holes, etc.) 

5.4.1. Matter response 

Cosmic structure, which growth in the matter era is driven 
by potential <£>, is very sensitive to modifications of <£> by new 
physics at low redshifts. As noted in Sec. 15.21 the Hubble fric- 
tion diminishes its sensitivity toward higher redshifts. The mat- 
ter power spectrum and other characteristics of cosmic structure 
are therefore the natural probes of the scenarios in the above 
classes (b) and (c). 

An example is shown on the right panel of Fig. [2] For both 
scenarios displayed on this panel the background is Einstein-de 
Sitter. The solid curves show the growing matter overdensity 
and gravitational potential in a pure CDM phase. [The plotted 
potential is weighted by t according to eq. d24l>.l The dashed 
curves display the same quantities for a toy scenario in which 
CDM constitutes only 50% of the energy density, while the 
other 50% is in a classical scalar field with w<f, = wcdm = 0. 
Since for the field Sp^/Sp^ = c 2 eS = 1, its clustering is sup- 
pressed and its perturbations do not contribute to the subhorizon 
potential. Correspondingly, the structure growth in the second 
scenario is seen to be significantly suppressed. Note that both 
scenarios have identical background expansion and the standard 
laws of gravity. 

Quantitatively, the impact of the potential $ over time (r,r + 
At) on matter overdensity d c observed at a later time r' ^S> r 
scales roughly as (tAt)$. The prefactor tAt quantifies the 
stronger suppression of an early impact by the higher Hubble 
friction, experienced by matter particles during the faster cos- 
mological expansion. The derivation follows straightforwardly 
from eq. d7): A nonzero potential $(r) over the time interval 
(r, r + At) contributes to the matter overdensity at a later time 
r'afl 

A(*£i) =Ad c {T') = -k 2 a(TmT)AT l ' lT ' 



9c ) r 



a(T") 



-r$(T)rAT. 



(24) 



In the models with standard gravity, CDM, and insignifi- 
cant clustering of dark energy, the dark energy perturbations 

" The last estimate in eq. {24} assumes that w < 1/3, so that o(r) grows as 
r or faster. The zero of conformal time r is chosen, as usual, as r — » when 
a — > 0. Logarithmic corrections that appear for w =1/3 (radiation domination) 
are ignored. 
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FIG. 2. — Growth of CDM overdensity d c (r,k) (rising black) during radiation dominatio n (le ft) and matter domination (right). Brown curves show the value of 
t$, responsible for the ultimate growth of Sp c /p c by a future time of its observation, eq. (24}. — Left: coupled photons only (dashed), photons plus 3 standard 
neutrinos (wide solid), and photons plus a tracking (w^ = 1/3) scalar field <j> replacing the neutrino density of the previous model (thin solid). — Right: pressureless 
matter only (solid), and a model with equal densities of matter and a scalar field that tracks it, = (dashed). 



can leave their mark on the potential only during horizon en- 
try. By eq. (l24l >. the corresponding early potential contributes 
to structure growth much less than the potential generated long 
after the entry by the clustered CDM and baryons. Thus in 
the absence of the non-standard physics of the types (b) or (c) 
above, it is justified to consider an approximate consistency re- 
lation between the backgr ound equation of state w(z) and the 
growth of cosmic structur e dKnox et al l2005Hlshak etal.l2 006: 
IChiba & Takahashill2007l) . 

As we see next, the situation is essentially opposite when 
the dark sectors are probed by the primary anisotropies of the 
CMB. 

5.4.2. CMB response 

Metric perturbations have a dramatic impact on CMB tem- 
perature anisotropies. In sharp contrast to large scale structure, 
primary CMB anisotropies (generated by linear evolution) de- 
pend most strongly on the values and evolution of gravitational 
potentials during horizon entry. They are only mildly sensitive 
to the potentials on subhorizon scales. (This need not to apply 
to secondary anisotropies.) 

The scalar perturbations of the CMB respond primarily to 
the sum of the Newtonian-gauge potentials $ + on all scales 
and all times (Sec. 13. It . While the inertial drag of the CMB 
by baryons, affected by $ alone, is important for the CMB sen- 
sitivit y to the baryon density and is noticeable around decou- 
pling dHu & Su givama 1996), it never dominates the evolution 
of CMB perturbations. For developing a general understanding 
of the CMB sensitivity to dark dynamics we will often ignore 
it. 

We stressed in Sec. [2] that dark inhomogeneities at a (con- 
formal) time t and the metric perturbations generated by them 
affect only the modes with k>H~ 1/r. The corresponding 
dynamics of CMB perturbations differs qualitatively in a tightly 
coupled regime, prior to hydrogen recombination at z lec ~ 1 1 00, 
and in a streaming regime, after recombination. We will con- 
sider the CMB response to the metric perturbation $ + "J 7 at 
these two epochs in turns. In the subsequent Sec. l5.4.3"l we will 
quantify this response by simple analytical formulas. We will 
see that despite numerous differences in the evolution of CMB 
modes before and after recombination, the CMB response to 



D = -3(0 + ¥ + R b O) 

| -3(0 + <rW (i) = kc s 

i Hubble and Silk damping 



i 30, = <L+3(<D + , F) 

eff ' 

FIG. 3. — An accurate mechanical analogy for general-relativistic evolu- 
tion (2) of an acoustic CMB mode. The CMB overdensity dj equals the de- 
noted distance to the tip of a pendulum with an internal frequency u) = kc s and 
its suspension point driven as specified by the gravitational driving term D, 
eq. 0. For adiabatic initial conditions the evolution starts with d m = 3£j„, 
where Q n is the superhorizon value of the Bardeen curvature. In the radiation 
era, initially, D m - 1(1 + 7 Rv)/Q- + TeRv )^jg an d m the matter era D m = | d m 
i Bashinskv & Seliak 2004), c.f. Figs.[T] a nd|3] The advantages of this approach 
over considering ST^^/T or eff = <5r (Ncwt) /r + $ = |d 7 as inde- 

pendent dynamical variables include: (A) Independence of the gravitational 
force that drives t/ 7 from d~, itself. (B) Epoch- and scale-independence of the 
relation between the superhorizon perturbation d m and the inflation-generated 
conserved curvature Q n (d[ n = 3Q n ). (C) Direct cause-effect connection be- 
tween local physical dynamics and the apparent changes of the variable that 
describes CMB perturbations. 



$ + \P at the respective epochs is very similar. 



Coupled regime fz > z rec ~ 1 100 J 

Any CMB mode prior to recombination evolves according to 
eq. (O of a driven damped harmonic oscillator. Fig.[3]presents 
an equivalent mechanical system — a pendulum whose evolu- 
tion is described by the same equation. The pendulum's internal 
frequency is kc s and its pivot is moved by a distance 

D=-3($ + V+R b $). (25) 

The CMB overdensity d 7 (r) then numerically equals the dis- 
tance to the pendulum's tip, as marked on the figure. To ac- 
count for the Hubble and Silk damping, the pendulum may 
be imagined submerged in a viscous fluid with appropriate 
time-dependent viscosity. All the influence of the metric per- 
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FIG. 4. — Gravitational suppression of CMB temperature anisotropy by growing cosmic structure in the matter era (5-fold for AT/T, 25-fold for power C/). — 
Left top: Evolution of d 1 and the driving term D during matter domination with R/, set negligible (solid). Evolution of d 1 from the same primordial perturbation if the 
metric becomes homogeneous before the horizon entry (dashed). — Right: Suppression of the CMB temperature power C/ for £ < 100 by CDM inhomogeneities. 
The solid curve shows C; in the concordance ACDM model with adiabatic initial conditions. The dashed curve describes the same model with changed initial 
CDM perturbations: cfcDM is artificially set to zero on superhorizon scales (the superhorizon d 1 , d y , and dt are unchanged). Then CDM inhomogeneities and the 
associated potential in the matter era are reduced. The smoother metric suppresses the CMB power at £ < 100 less, as the C; plots (obtained with CMBFAST) 
clearly show. — Left bottom: The plots of 4> and \t transfer functions confirm that the last model has smaller <E> + \t, hence, the enhancement of its power at low £ 
cannot be attributed to the ISW effect. 



turbations on d~ 1 is transferred through the (generally, time- 
dependent) position of the pendulum's pivot D(t). 

This mechanical analogy of the acoustic CMB evolution 
makes intuitive the following important facts: Foremost, the 
CMB overdensity is influenced strongly by the values of the po- 
tentials at horizon entry. Also, as already observed in Sec. 15.21 
the sensitivity and even the sign of the CMB response to the 
subsequent changes of the potentials depend strongly on the 
time and duration of the change. 

When at horizon entry the driving term D is close to the su- 
perhorizon value of dry (as it usually is for adiabatic pertur- 
bations) and when D does not decay quickly during the en- 
try then CMB temperature anisotropy is suppressed dramati- 
cally. For a specific example, roughly reflecting the evolution 
between radiation-matter equality and decoupling (3 x 10 3 > 
z > 1 x 10 3 ), let us ignore the baryon-photon ratio Rt, by tak- 
ing D ~ -3(<I> + *). Let us evaluate the potentials in a CDM- 
dominated limit, in which ^(r) = W(r) = const = — hd^, where 
for adiabatic perturbations d m = d a (T—>0) = 3((t— >0) for all 
species a[3 In these approximations, 

D(t) = -<4, (26) 

being time-independent during and after the mode entry. Then 
the amplitude of a coustic CMB oscillations is suppressed 5-fold 
dBashinskyil2006T) . which is evident from Fig. 0or the left top 
panel of Fig. 2] 

After the entry, a slow variation of D over a characteristic 
time St that exceeds the period of d 1 oscillations (Srk 3> 1) 
has minor impact on the oscillation amplitude and phase. In 
particular, this impact vanishes in the adiabatic limit, StIc^oo. 
A typical temporal scale of linear variations of <£> and "f, and 
so of D, is St ~ H~ l . This applies even to subhorizon scales, 

12 The CDM- domination result <E> = \& = — ^ d m follows straightforwardly from 
eqs. f7l ll U12l and from $ = ^ by eq. 4 1 4 1 . It is also easy to derive from the 
known relation between potentials and the traditional proper overdensity by 
remembering that for pressureless matter d = Sp/p—3^. 



where only the perturbations with negligible pressure, hence no 
rapid internal oscillations, contribute to the potentials. Thus 
deeply subhorizon CMB modes (for which Srk ~ k/Tt 3> 1) 
are insensitive to such variations. Conversely, the modes that 
are closer to horizon entry are affected more strongly. 

Streaming regime (z < z rec ~ 1 100) 

Analogous results exist for the sensitivity of the CMB to gravi- 
tational potentials after recombination. 

When the universe is matter-dominated and baryons have de- 
coupled from the CMB (10 3 > z > 1) then $ and * are equal 
and time-independent on all scales that are sufficiently large to 
evolve linearly and to be unaffected by the residual baryonic 
pressure. For time-independent $ + ^ and negligible collision 
term Cj in eq. ©, an "effective" CMB intensity perturbation 

teff = 6+4($ + 10 (27) 

is constant along the line o f sight (Kodama & Sasaki 1986; 
iHu & Sug iyama 19 94lll995h . Under a general evolution of the 
potentials, eq. ((5]) gives 

4ff+n,V,i e ff = 4($+*) + C r . (28) 

Prior to horizon entry, t he intensity pert urbation i is frozen 
and equals i m = |t/ 7 ,i n (see lBashinsky||2006l) . As discussed ear- 
lier, for the modes that enter during the domination of pressure- 
less matter, $(r) = \&(t) = — i<4i- Then, for adiabatic perturba- 
tions for which c/ 7 ,j n = d m , eq. (f27l > gives 

ieff = -^in- (29) 

Similarly to the CMB evolution in the coupled regime, for 
most of these modes the late ISW effect due to the slow de- 
cay of the potentials at z < 1 does not noticeably change i e ffH 
Thus, first, the presently observed temperature anisotropy in the 

13 The solution of the transport equation )28t for a perturbation mode t c ff = 
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considered modes is suppressed by the potentials of the clus- 
tering matter 5-folc0 — suppressed with respect to a fictitious 
scenario in which in the matter era $ = ^ = 0, consequently, 
teff = '-in- Second, after decoupling as well, the primary CMB 
anisotropics have little sensitivity to the potentials' evolution on 
subhorizon scales. 

The 5-fold suppression is a physical effect, independent of 
the choice of variables or gauge. It is absent in models where 
the matter-era potentials decay early during horizon entry or 
where the laws of gravity are modified so that matter inhomo- 
geneities do not perturb the metric. Note that if all modes that 
contribute to the CMB temperature autocorrelation Q at some 
£ are suppressed 5-fold then the Q is suppressed by 5 2 = 25 
times. 

The right panel of Fig. |4] presents convincing evidence for 
the reality of the order-of-magnitude suppression of the CMB 
temperature power spectrum at low fs by dark matter cluster- 
ing. The plot shows the CMB temperatu re power spectra C/ 
obtai ned with a modified CMBFAST code (ISeliak & Z aldarriaga 
1996*1) for two models both of which have the identical matter 
content of the concordance ACDM mo del dSpereel et al.l l2006; 
Selia k et alj20 06: Teg mark et al.l2 006). The models differ only 
by the initial conditions for CDM perturbations. In one of the 
models (solid curve) all species are initially perturbed adiabat- 
ically. In the other model (dashed curve) the CDM density 
perturbation d c is artificially set to zero on superhorizon scales 
while the initial values for d l7 d v , and dt, are unchanged. As a 
consequence, in the second model the metric inhomogeneities 
in the matter era are reduced. While this model has a smaller 
ISW effect (see the left bottom panel of Fig. 3), due to the 
smoother metric, its CMB power for I < 100 is considerably 
larger. 

The suppression of the CMB anisotropy on large scales, 
which enter during matter domination and active growth of 
structure, should not be partly traded for the "resonant self- 
gravitational driving" of small-scale modes that enter in the ra- 
diation era. The acoustic modes entering during radiation domi- 
nation are often said to be resonantly driven by a specially timed 
decay of their self-generated potential. In reality, the suppres- 
sion of large relative to small scales is not very sens itive to the 
timin g of the potential decay in the radiation era (Bashinsky 
2006). Moreover, the same (when accounting for neutrinos, 
even somewhat larger) suppression of large relatively to small 
scales would be observed if in the radiation era the metric were 
unperturbed, hence, the sm all-scale modes o bjectively could 
not be driven gravitationally (Bashinsky 2003). 

5.4.3. Quantifying the response of the CMB power spectrum 



A(t) 



,zk-(x-nT) : 



in a potential $ + ^ = F(t) e' k * is 



A(r)=A in + 



4 [ T F(r')e ik nT 'dT'. 
Jo 



(30) 



Due to the oscillating factor e' k nT , any slow variation of the subhorizon po- 
tential over a temporal scale St does not affect any modes except for those few 
with |k ■ n| St < 1, i.e., those whose wavevector is nearly orthogonal to the 
direction of light propagation. 

14 This suppression, by a factor of 5, is larger than an apparent factor of 2 sup- 
pression that is usually inferred, erroneously, from the traditional description 
of the Sachs-Wolfe effect in terms of the proper Newtonian-gauge perturba- 
tions. (Indeed, for the latter AT/T| observed ~ -i(5r< Newt »/7-) d e C0up i i ng , assum- 
ing matter domination at decoupling. Yet this result is specific to the Newtonian 
gauge condition for the decoupling surface, on which the considered perturba- 
tions are superhoriz on. It should not be given its traditional physical interpre- 
tation directly.) See Bashinsky (2006) for additional discussions. 



The CMB sensitivity to gravitational potentials at var- 
ious epochs can be quantified as follows. The tempera- 
ture anisotropy of CMB radiation observed in a direction 
n is given by the followin g integral along the line of sight 
dSeljak & Zaldarriagall996l) : 

- 1 dr S(r,nr(r)). (31) 



Here, To is time at present, r{f) = tq-t (in flat models) is the 
radial distance along the line of sight, and the source equals 



S = g[^d 7 + $ + f - vim + Q'n-flj I + g I $ + * 



(32) 



with v' b being the baryon velocity, Q'j being determined by the 
radiation quadrupole and polarization, and a visibility function 
g(r) = exp(- dr/Tr) giving the probability of CMB photons 
to reach us unscattered. Regardless of the random realization of 
the primordial curvature perturbation Q n , the physics behind the 
subsequently developed perturbations of potentials, densities, 
etc. in linear theory can be described by transfer functions 

7*(fc,T) = $(fc,T)/CirA'*), (33) 

with analogous definitions for other scalar perturbations: 'J?, d 1 , 
velocity potential ut, (s.t. v' b = -diut,), etc. 

The CMB constraints on the transfer functions are derived 
from observational estimates of CMB angular power spec- 
tra C/. In particular, the te mperature power spectrum equals 
dSeliak & Zaldarriagall996h 

2 



dk 

T 



A 2 (k) 



drT s (r,k)ji(kr(T)) 



(34) 



k 3 P((k) / (2ir 2 ) is the dimensionless primordial 



Q = 4tt 
where A 2 c (k) : 

power of CinW, an d ji is the spherical Bessel function. In 
eq. (1341 . the full transfer function T$ of the source S, eq. d32l . 
is a differential operator, with derivatives corresponding to the 
direction-dependent terms of S. In the following discussion we 
will mostly be concerned with its scalar terms, as only they in- 
volve the gravitational potentials directly. 

It is useful to keep in mind that the dominant contribution to 
Ci at a given i comes from the modes with k ~ l/r. Indeed, 
ji(kr) vanishes exponentially at smaller values of k and as 1/k 
at higher values. 

Radiation era (z > z eq ~ 3000, probed by l> 200) 

By the arguments of Sec.|2j only the modes that enter the hori- 
zon before radiation-matter equality provide information about 
the metric and dark dynamics in the radiation era. Then af- 
ter horizon entry dp/ p w 8p 1 / p 7 oscillates with nearly con- 
stant amplitude. Consequently, the induced $ + 'J? decays as 
r" 2 oc a~ 2 [eq. (fl)]. Although the potential decay is eventu- 
ally halted by the growth of structure since matter domination, 
$ + 'J? for these modes remains a subdominant source of CMB 
anisotropy in eq. 1321 . exceeded by the intrinsic photon-baryon 
perturbationsfij 

Thus the impact of gravitational potentials in the radiation 
era is confined to horizon entry. It is fully encoded in the ampli- 
tude and phase of the subsequent acoustic oscillations. Specific 
cases of such an impact were considered in Secs. l5.ll and l5.2l 

15 Photon- baryon perturbations al s o decay on sma ll scales due to Silk damping. 
Estimates i Hu & Sugivama 1996; Weinberg 2002) show that the exponential 
Silk damping overcomes the quadratic decay of the potential only at I > 4000, 
where secondary anisotropies are expected to dominate the considered primary 
signal. 
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From equality to recombination (z rec < z < Z eq , 100 < i < 200) 

The CMB modes that enter during and after radiation-matter 
equality (z eq ~ 3000) oscillate in a significant gravitational po- 
tential of growing matter inhomogeneities. Close to z e q, the 
driving potential <£> + '!' still evolves appreciably. The poten- 
tial continues to decay because of the residual radiation den- 
sity. It decays until recombination (z rec ss 1 100) also due to cou- 
pling of the baryonic component of matter to the CMB, slow- 
ing the structure growt h by excluding the baryons from it (e.g., 
iHu & Sugiva"malll996l) . 

Altogether, these effects lead to complicated evolution of the 
modes that enter during z rec < z < z eq - Then both the intrinsic 



and gravitational terms in the source of AT /T (l32l contribute 
noticeably. As highlighted by many studies, the term $ + ^ in 
eq. d32l due to the continuing decay of potentials then is also 
large and boosts the height of the first acoustic peak (early ISW 
effect). 

After recombination (z < Zrec — 1 100, probed by £ < 100) 

The evolution of linear perturbations becomes simple again 
when z <C z re c- On the scales that enter at this epoch but before 
dark energy becomes dynamically relevant, we expect that dur- 
ing linear evolution $(r) = $(r) = -^c/; n = const (footnote[T2b. 
The presently observed perturbation of CMB intensity for al- 
most all these modes is suppressed 5-fold. Our next goals will 
be, first, to confirm that the contribution of these modes to Q is 
indeed gravitationally reduced by a factor of 5 2 = 25. Second, 
to find the sensitivity of Q to different values of potentials due 
to non-standard physics. And third, to quantify the Q response 
to the evolution of potentials on subhorizon scales. 

We start from noting that not only do the modes of differ- 
ent k's contribute to the correlation function Cg ( 134b incoher- 
ently, but coherence is also lost for same-fc contributions at 
widely separated times, |r — r'\k ^> 1. To see this explicitly, 
we rewrite eq. d34l i as 

Cp = A^J^A 2 {k) JdrT s (T,k) JdT%*(T',k) p £ (kr,kr')- (35) 

The last factor pe(kr(T),kr(T')) describes the kernel of project- 
ing the harmonic plane-wave modes on a spherical multipole I 

pi(x,x') = ji(x)jt(x'). (36) 

This function for I = 10, as an example, is shown in Fig. [5] 

In the integrand of eq. d35l ), except for very low redshifts, 
Pi(kr,kr') varies with k much more rapidly than Ts(r,k) and 
r/(T',/fc)|3 With T s , T s *, and A 2 ( (k) almost unchanged over 
many k periods of the pi oscillations, we can see from Fig. [5] 
that positive and negative contributions to f(dk/k) pi(kr,kr') 
mutually cancel whenever 

\t-t'\ = \r-r'\ > 1 /(contributing k). (37) 

Thus the contributions to C/ from sources at this temporal sepa- 
ration are incoherent. In particular, we can ignore the coherence 
of, and study independently, the contributions to C; at the hori- 
zon entry (r < 1 jk) and during subhorizon evolution (t> l/k). 

First, we consider a time interval from r = to r ent ~ l/k. 
For a mode that enters at any time after recombination, the 

16 The characteristic scales of variation are Ak ~ min(r~' , r'" 1 ) for pi, versus 
Ak ~ 7~ and t' _1 for 7$ and T? respectively. 
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FIG. 5. — The isocontours of the ker nel pi( x\ ,x 2 ) = jt(xi)jt(xi) of pro- 
jecting perturbation &-modes to Q, eqs. J35II361 . The figure is for i = 10. The 
labels on the contours show the values of pio/lCT 3 . Note that pt{x\ ,x 2 ) ~ 
whenever either x\ < I or Xz < I. Also note that pi>0 along the line x\ = x 2 
(diagonal dashed line), but pi oscillates through positive and negative values 
along any line X\ = cx 2 with c ^ 1 (e.g., lower dashed line). 



source (|32l forms a complete derivative Q 



d_ 



-t/ 7 , in + $ + * 



(38) 



Starting from eq. d34l ), neglecting the change of ji(kr(r)) over 
the considered time interval (0,r ent ), and trivially integrating 
the complete derivatives d38b over dr and dr', we obtain 

2 



sc. 



,(postrec entry)^ 



47T 



dk 



A 2 c (k)jf(kr) 



-d lM + <S>+^ 



where the variables d-, 



(39) 

entry 

$, and \P now stand for the corre- 



sponding transfer functions, normalized by the condition Q(t — > 
0, k) = 1, When the primordial power is nearly scale invariant, 
we can neglect the ^-dependence of A?(&) over the range of 



k ~ £/r that contribute to the integral d39l . We can also ignore 
the variation of the transfer functions in the last parentheses 
over this range. Integration of the remaining ^-dependent terms 



by the standard formula / — jf(x) = [2£(£+ l)]" 1 gives 



sc ( r 



strec entry)^ 



2wA 2 ( (£/r ellt ) 
£{£+!) 



-d 7(in + $ + * 



(40) 



entry 



Eqs. ( 1391 or (l40b quantify the response of C\ to the gravita- 
tional potentials at horizon entry after recombination. These 
equations confirm that for adiabatic perturbations that enter in 
the matter era the potentials $ = = -^d m suppress the corre- 
sponding contribution to C/ 25-fold. 

We now consider the response of C; to changes of potentials 
on subhorizon scales. Then the C MB sensitivity can be quanti- 
fied in the Limber approxima tion (Limber 1954), applied to the 
CMB bv lKaiserl dl984lll992l) and lHu & White! d 19961) . With the 
details of the calculation given in Appendix[A] the response of 
the power spectrum of CMB temperature to subhorizon changes 
of potentials is roughly 

2ir 2 A 2 (£/r) . r 

17 For such modes, g significantly deviates from zero only when the mode is 
superhorizon. Then the factor that multiplies g in the full source (32) reduces 
to (i(i ri i n +$ + *), giving eq. (38). 



j£-i(postrec subhor.) 



(41) 
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FIG. 6. — While the acoustic pattern of CMB temperature anisotropy is 
formed on the surface of last scattering (z rec ~ 1 100), the metric perturbations 
on this surface do not play a major role in the observed CMB spectra Q . The 
potentials at last scattering do control the subdominant "baryon loading" ef- 
fect. Primarily, however, the CMB multipoles are sensitive to the gravitational 
potentials and underlying dark dynamics at the horizon entry of the corre- 
sponding modes (z ~ £ 2 /4 for z. < Zeq, and z. ~ 20 ^ for z Zeq), covering the 
entire range from z ~ 1 for lowest I up to z ~ 6 X 10 4 for I ~ 3000. 



Here r is the comoving radial distance to the changing poten- 
tials and St is the time taken by a change <5(<i> + , i'), assumed 
gradual. 

Even for large changes of potentials (|£(<i> + \E r )| ~ 1$ + *!) 
the observed response 6Cf uhhol ) d4lT> is suppressed relative to 
C^ entry) d40b by a factor We argued previously that un- 
der linear evolution generally 5t ~ fi~ l . Then for subhorizon 
scales (r/l <C ri~ x ) the factor is much smaller than unity. 
This suppression is easily understood as the cancellation of pos- 
itive and negative contributions to AT /T from the peaks and 
troughs of all perturbation modes ex cept for those whos e k is 
almost orthogonal to the line of sight (lHu& Whitdll996h . 

In addition, dark dynamics influences the CMB indirectly 
through the gravitational impact on baryons, coupled to the 
CMB by Thompson scattering. After recombination, the 
baryon-mediated linear contribution on subhorizon scales is 
suppressed even more than the late ISW contribution ( f4~LT >. 
Indeed, the direct contribution of baryon velocity v l b n; to the 
source of CM B temperature per turbation d32l in the Limber 
limit vanishes dHu & White! 19961) : In this limit, only the modes 
with k orthogonal to n would contribute to Q, yet for their 
scalar perturbations (with v' parallel to k) v' b rij = 0. The indi- 
rect effect of weakly coupled baryonic perturbations through 
affecting the photon overdensity or other photon multipoles in 
the source ( f32b is also diminished relative to the already sup- 
pressed ISW term by an additional factor St/tt, where tj is the 
photon free-flight time. This factor is small for any temporal 
interval after CMB decoupling. 

The CMB is also affected through various non-linear effects, 
e.g., its lensing by cosmic structure or Sunyaev-Zeldovich ef- 
fect. Unlike the primary anisotropics, these non-linear features 
may be useful probes of the dark dynamics on small scales 
whenever they can be separated from the dominant primary sig- 
nal (e.g., Seliak & Zaldarriaga 1999; Hu & Okamoto 2002). 



tials during recombination are of relatively minor importance 
for the observed anisotropies. (Exceptions are the scales of the 
first peak that happened to enter during recombination, and the 
nondegenerate baryon-loading signature on smaller scales.) 

When the potentials do not decay quickly during the entry, 
they significantly suppress the CMB temperature spectrum Q 
(by a factor of 25 in the Einstein-de Sitter scenario with adia- 
batic perturbations.) The Q response to changes of the poten- 
tials after the entry is weak. After recombination, in particular, 
the Q response to the subhorizon changes of potentials is di- 
minished by a factor ~ H/k. 

In principle, AT /T of the CMB responds comparably well 
to the metric inhomogeneities during the horizon entry at either 
high or low redshifts. Nevertheless, the CMB is a consider- 
ably better probe of the horizon-scale potentials at high red- 
shifts because of the reduced cosmic variance. Let us define the 
wavevector k ent of the CMB modes that enter at a redshift z by 
a condition k ent S{z) = 1, where S = J c s cIt is the corresponding 
size of acoustic horizon. The modes that enter at a redshift z 
match to the multipole £ ent = k eM r = r/S. Taking S ~ t/vj as a 
valid estimate both in the radiation era (when R/, is negligible) 
and after recombination (when baryons decouple and formally 
Cg ff = 1 /3), we obtain 



*(z) : 



z/18, radiation era, 
~ 2^/z, matter era, 



(42) 



estimated for a ACDM model with Q m h 2 = 0.13 and h = 0.7. 
For this model the redshift of equality z eq « 3 1 10 and recombi- 
nation z rec ~ 1090 match to £ ent ps 200 and 90 respectively. 

The 1-er uncertainty of probing a non-degene rate effect, 

quant ified by a parameter p, can be evaluated as (e.g jKendall & Stuartl 
fl979l) 



Ap= 1 



v ( dQ/dp 



(43) 



For a CMB experiment limited on the studied scales only by 
cosmic variance, for either temperature o r polarization power 
spectra (r.m.s. C,)/C, = \(£+ ±)f sky r 1/2 (iKnoxl 11995b ISeliakJ 
1996l lZaldarriaga et alj 19971) . where / s k y < 1 is the experimen- 
tal sky coverage. Then 



d\nQ 
~dp~ 



(44) 



Thus for constraining an effect that affects the modes which 
enter the horizon at a redshift z over a span of redshifts Az, by 
eq. 1421 . we can expect accuracy 



Ap 



•f^idhyQ/dpT 



radiation era, 



matter era. 



(45) 



This accuracy improves with increased Az and, when probing 
the radiation era, with increased z. 



Overall 

The above results can be summarized concisely as follows. At 
a given comoving scale k, probed by I ~ kr w k x 14Gpc, the 
CMB anisotropy is most sensitive to the value and evolution of 
$ + ^ during the scale's horizon entry. This is illustrated by 
Fig-H] Although most of the observed CMB photons scattered 
last during hydrogen recombination at z rec ~ 1100, the poten- 



6. APPLICATIONS AND CAVEATS FOR MODIFIED GRAVITY 

We now briefly consider the possibility that general relativ- 
ity (GR) breaks down on cosmological scales. Modified gravity 
(MG) offers rich phenomenology by invoking new degrees of 
freedom whose dynamics is substantially different and often in- 
volves more parameters than the standard cosmological model. 
It typically predicts nonstandard evolution of perturbations on 
horizon and subhorizon scales alike. 
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While the detection of such phenomena as non-standard 
growth of cosmic structure or anomalous lensing may indicate 
MG, we will see that without any restrictions on the dark dy- 
namics, the identical effects could always be generated by a 
non-minimal dark sector that influences the visible matter ac- 
cording to the standard Einstein equations^ Further in this 
section we will argue that other features should yet allow to 
discriminate MG observationally. 

We will assume that even if full Einstein gravity fails on cos- 
mological scales, the Einstein principle of equivalence remains 
valid for the visible species. This assumption is common to 
many existing MG models. It is motivated by the relatively 
strong terrestrial and solar-system constraints on the equiva- 
lence principle. 

Thus we suppose that the regular matter couples covariantly 
to a certain matter-frame "physical" metric g^ u . However, we 
now neither take for granted that all dark fields also couple co- 
variantly to the same metric g M „, nor assume that the dynamics 
of g^ u itself is governed by the Einstein equations. 

Under the weaker assumption of the equivalence principle 
for only the visible matter, all observable signatures of new 
physics can still be quantified by any of the three parameter- 
ization schemes of Sec. [4] Indeed, the g M „ background can still 
be described by a single number for its present spatial curvature 
and by its uniform redshift-dependent expansion rate H (z). The 
potentials $ and \I>, defined by eq. (O to parameterize the in- 
homogeneities of the physical metric will play their usual 
role in the evolution of light and baryons. Moreover, the (effec- 
tive dark) energy and momentum densities assigned to the miss- 
ing sources of curvature by the naive application of the Einstein 
equations will evolve in agreement with the usual local conser- 
vation laws, which is easily seen as follows. 

Let by definition 

t^=^g^- y j rr, (46) 



8ttG 



known a 



where the last sum is over the known regular particles. Their 
energy-momentum tensor [constructed unambiguously from 
the species' action S a as T£ v = (2/ v / = g) SSa/Sg^] is covari- 
antly conserved by the assumed covariance for the regular 
species. The matter-frame Einstein tensor G^ v = R^ v - ^g^R, 
where is the Ricci tensor of the physical metric, is also 
covariantly conserved by the Bianchi identities. Thus the entire 
expression (l46l l is covariantly conserved: 



1 eS &drk ;u 



(47) 



with all covariant derivatives being taken relative to the physical 
metric. 

Since T^ dmk is covariantly conserved, the background and 
perturbations of the missing energy and momentum evolve ac- 
cording to equations ( fToT l - ( [T9| >, derived from the identical con- 
servation law. Thus if all our probes of the invisible degrees of 
freedom are based solely on their gravitational impact on light, 
baryons, and other regular particles (neutrinos, WIMP's when 
probing dark energy, etc.) then phenomenologically all observ- 
able signatures of a considered MG model can be mimicked 
by an effective GR-coupled dark sector. Specifically, we can 
find the corresponding effective w(z) to reproduce the missing 
energy background, and the effective anisotropic stress a(z, k) 
and stiffness cl s (z,k) to describe scalar perturbations. For ex- 

18 We stress that these features remain useful hints of MG. Since they are even 
more ubiquitous than the scenarios of modified gravity and may reveal other 
new physics, they are well worth searching for. 



ample, for a popular DGP gravi ty model (Dvali et al.l2 000) this 
was demonstrated explicitly bv lKunz & SaponeT i 20071) . 

With this discouraging general conclusion, we may enquire 
whether cosmology can at all reveal definitive distinctions be- 
tween MG and general relativity with a peculiar yet physically 
permissible dark sector. To establish such distinctions, let us 
summarize the conceptual differences of MG from general rel- 
ativity: 

A. Some dark degrees of freedom may not couple covari- 
antly to the matter metric g^ v . 

B . The gravitational action may not be given entirely by the 
Hilbert-Einstein term S grdv = (16ttG)~ 1 J d x^J^R. 

The distinctive observable consequences of these special prop- 
erties of MG may include: 

1. The effective dark dynamics, which is observationally 
inferred by assuming the Einstein equations, violates 
the equivalence principle (EP). The EP violation can be 
seen, e.g., as 

i. The dependence of the inferred local dark dynam- 
ics on the distribution of visible matter, when it 
cannot be explained by non-gravitational dark- 
visible coupling allowed by particle experiments. 

ii. Superluminality of the inferred dark dynamics. 

2. The dynamics of gravitational waves (tensor modes) de- 
viates from the predictions of the Einstein equations, as- 
suming that both the visible and inferred dark species 
contribute to the energy-momentum tensor in the sim- 
plest way. 

Most of these signatures have already been utilized for fal- 
sifying MG mode l s with existing or sugg ested observations, 
(Clo we et alj2006HBradac et al.ll2006l) and dKahva & Woodardl 
120071) : we comment additionally on them next. 

6.1. EP violation for the inferred dark dynamics 

The violation of the first condition can be illustrated by an 
extreme toy theory in which the regular matter ("baryons," for 
short) constitute all independent degrees of freedom. Let the 
metric in this theory be specified by baryon distribution via 
some deterministic relation (e.g., as g^ = A _I r b ^ on , follow- 
ing from an action S = Sbaryon ~ §d A x^ r gk with A being a 
constant). Even in such a contrived theory, by the above ar- 
guments, the effective missing energy and momentum densi- 
ties d46b would appear to evolve and gravitate in agreement with 
energy-momentum conservation and the Einstein equations. In 
this example, however, the effective dark density and stress are 
uniquely determined by the distribution of the visible matter. 
This does not occur for truly independent dark degrees of free- 
dom. 

In more realistic MG theories we should not expect a deter- 
ministic relation between the visible and effective dark distri- 
butions. Still, if the dark and visible sectors interact other than 
by coupling covariantly to the common metric then the inferred 
laws of the effective local dark dynamics would depend on the 
visible environment. 

Detection of such dependencies would be particularly feasi- 
ble for the dark matter, for which there are plentiful observ- 
able regions with varying environment: varying in both visi- 
ble matter density and in its ratio to dark density. In addition 
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to baryons clumping strongly at low redshifts, the density of 
all known visible species varies with redshift due to the Hub- 
ble expansion. The ratio of visible to dark density is perturbed 
even on large scales due to dark matter being decoupled from 
the photon-baryon acoustic oscillations and clustering on its 
own until recombination. The segregation of dark matter and 
baryons is even more apparent at late times on the scales of clus- 
ters and smaller, be coming almost complete in the famous bul- 
let cluster example dMarkevitch et alJl2004t IClowe et~alT 2006: 
iBradac et alj|2006l) . 

In order to test that the standard CDM+GR model accounts 
for the signatures of dark matter at various redshifts and scales, 
it is crucial to observe and robustly model the dynamics on a 
wide range of scales, including the highly non-linear ones. It is 
also important to utilize complimentary probes, such as bary- 
onic matter (responding to $), and the CMB or gravitational 
lensing (probing $ + Various signatures of dark matter that 
allow to compare dark matter parameters under different condi- 
tions include: the height alteration of the odd and even CMB 
acoustic peaks due to the CDM potential affecting the cou- 
pled baryons before recombination dHu & Sugivamal 19961). th e 
early ISW enhancement of the first peak (e.g.. lHu et al.ll 19971) . 
the significant suppression of the CMB temperature power C; 
below the first peak (Sec. 15.4.31 and Fig. |4), the linear and 
nonlinear dynamics of cosmic structure, and the lensing of the 
CMB and background galaxies by the CDM potential at lower 
redshifts. 

Dark energy (DE) does not trace any visible species at low 
redshifts, when its density appears almost redshift-independent. 
DE may track the radiation or matter background at higher red- 
shifts; then its perturbations are not expected to evolve similarly 
to any of the standard species. (For example, see the evolution 
of the perturbations of tracking quintessence in the radiation 
and matter eras on the right panels of Figs. Q] and [2] respec- 
tively.) Constraints on the background equation of state w can 
rule out specific DE or MG models but do not decisively dif- 
ferentiate between the DE and MG paradigms. Beyond w, the 
perturbations of dark energy are likely to be stiff and thus can 
be constrained only during horizon entry. As discussed previ- 
ously, the CMB can provide such constraints over a wide range 
of redshifts, from z ~ 1 up to z ~ 10 5 , becoming much tighter 
toward higher redshifts (Sec. [5] and Fig.|6j. 

With the cosmological constant fitting the current data well, 
the search for any manifestation of non-trivial DE or MG dy- 
namics may be the utmost priority for establishing the origin of 
cosmic acceleration. This includes falsifying the background 
equation of state w(z) = — 1 at low redshifts but by no means 
is limited to w(z). Examples of other searches, which may 
end up being more fruitful, are: the search for a subdomi- 
nant non-standard component in the radiation or early matter 
era ( Sees. 15.21 and 15.41 ), nonstandard grow th of cosmic struc - 
ture (Linder 2005), and nonstandard dBertschingerll2006l) . 
probed by comparing the CMB or lensing signal with the galaxy 
and cluster distributions. 

6.2. Superluminal dark flows 

Another possible signature of MG is the superluminal prop- 
agation of dark perturbations. The superluminality may be spu- 
rious, from our describing MG by a parameterization that as- 
sumes general relativity. It may also be real, from the equiva- 
lence principle not applying to some dark degrees of freedom. 

While numerous observational tests for superluminality can 
be thought of, here we consider only one specific example. Spu- 



rious or real superluminal propagation of (effective) dark inho- 
mogeneities can be constrained by the sensitivity of the CMB 
acoustic phase to the propagation velocity of dark perturbations 
(Sec. 15.31 ). A straightforward Fisher-matrix forecast shows that 
TT, TE, and especially EE spectra from the future high-^ CMB 
experiments, Planclo and AClFl in particular, will strongly re- 
strict the abundance of any dark component that supports scalar 
perturbations with c 2 p > 1/3, including c p > 1. 

If fi-polariza tion of the CMB reveals the signal of relic gravi- 
tation al waves dSeliak & Zaldarria gall997l : lkamionkowski et al.l 
1 1997b . a similar effect in the tensor sector will directly con- 
strain the streaming dark species for which streaming velocity 
c p exceeds the speed of gravitational waves, c ! grav . [Indeed, 
c P . v > Cj.grav, albeit under so mewhat different condition s, for 
the dark matter emulator of Kahva & Woodardl d2007h . dis- 
cussed below.] The oscillation amplitude of the gravitational 
(tensor) modes after h orizon entry is noticeably affected by 
neutrino perturbations (Weinberg 2004). However , the phase 
shift of the tensor oscillations is strictly forbidden (Bashinsky 
120051) when for all species c p c iaik < c Sigrav , as required by GR. 
If, on the other hand, the velocity of neutrinos or other dark 
species with non-negligible a bundance exceeds c s , grav then the 
arguments of Bashinsky (2005) show that the phase of the ten- 
sor BB signal in the CMB will be necessarily shifted. 

6.3. Nonstandard phenomenology of gravitational waves 

An interesting test of a broad class of MG alternatives 
to da rk matter was recently suggested by iKahva & Woo dard 
(2007). They considered the propagation of a fundamen- 
tal tensor field g^ [whose kinetics is governed by S gmv = 
( 1 6irG)~ l J d 4 x^/^gR] in any of the model where g^ v is sourced 
by the luminous matter alone by the Einstein law 

(48) 

The physical metric g^ v was set to reproduce the observed grav- 
itational potentials in the vicinity of our galaxy. The knowledge 
of the other s pecifics of the MG dynam ics was therefore unnec- 
essary. Then K ahva & Woodardl (12007) showed that the arrival 
of the gravitational waves from a cosmological event, e.g. a su- 
pernova, is noticeably delayed with respect to the arrival of the 
associated neutrinos and photons. 

Wh ile the assumption d48l i. implicit in IKahva & Woodardl 
d2007l) . incorporates many models it is not generic. For ex- 
ample, it does not apply if 



^ 7r( ^ ^luminous' 



(49) 



where G differs from the local value of Newton's gra vitational 
constant G. Even in the TeVeS model (Bekenstein 2004), moti- 
vating (l48l l. G^G strictly, althou gh it is natural to a ssume that 
for TeVeS this difference is small (Bekenstein 2004). Yet more 
generally, G may depend on the new MG degrees of freedom, 
e.g., on the scalar field cj> of tensor-scalar or tensor- vector-scalar 
models. 

In any case, even if eq. (l48l i f ails and the exact pre diction for 
the gravity wave delay by Kahya & Woodard (2007) does not 
apply, the gravitational waves and neutrinos or photons would 
still generally propagate with different velocities. Thus the fu- 
ture gravitational wave astronomy can offer robust tests for dis- 
criminating GR-coupled dark matter from modified gravity. 



" http://www.rssd.esa.int/Planck 

20 http://www.physics.princeton.edu/act 



Sergei Bashinsky 



17 



7. SUMMARY 
7.1. Approach 

Constraints on the inhomogeneous dynamics of the dark sec- 
tors are essential to full understanding of their nature. Varying 
with both redshift and spatial scale k, the observable imprints of 
dark perturbations provide plentiful information about the dark 
sectors' local kinetic properties and (self) interactions. This 
information significantly complements the dark species' back- 
ground equation of state w(z). 

Reliable extraction of the information encoded in the dark 
perturbations is hindered by the indirectness of their, some- 
times subtle, gravitational impact on observables. It is also 
obstructed by the numerous contributions to the observables 
from other complex multiscale cosmological and astrophysi- 
cal phenomena. Moreover, controlled repeatable measurements 
of dark sectors cannot be afforded when the measuring device 
is the entire observable large-scale universe. Nevertheless, the 
unavoidable "nuisance" phenomena and contaminations can be 
counteracted by detailed understanding of the involved physics 
and by complementary probes of as many independent charac- 
teristics of the dark sectors as can be observationally accessed. 

For a tractable and systematic study of the dark sectors be- 
yond the background equation of state, we map the local kinetic 
properties of inhomogeneous dark dynamics at a given redshift 
to the characteristics of the observed cosmological distributions 
that can reflect those dark properties. The correct mapping is 
more likely to be achieved with a description of evolution that 
is simple and that manifests the objective causal relations ex- 
plicitly. 

The dynamics of gravitationally coupled perturbations of 
dark and visible species during horizon entry is pivotal to prob- 
ing the dark inhomogeneities. Then any dynamical species, in- 
cluding dark radiation and all dark energy candidates with w ^ 
-1, are necessarily perturbed and leave imprints of their highly 
specific inhomogeneous kinetics on the observable probes. The 
contribution of subhorizon perturbations of a given magni- 
tude Sp/p to gravitational potentials is suppressed, as $ ~ 
(H/k)Hp/p. 

On the scales of the horizon and beyond, the apparent gravi- 
tational impact of the dark species on the visible ones depends 
strongly on the description used. In most descriptions, the ap- 
parent inhomogeneities of the visible species are changed by 
the properties that dark species have long before and even long 
after the change. This misguides the identification of the ob- 
servable features that reflect the internal dark properties at spe- 
cific epochs. 

The apparent cause-effect mismatch is, nevertheless, not in- 
trinsic to linearly perturbed cosmological evolution. Within 
Einstein gravity, the internal local properties of dark species at 
a past time r affect only the perturbations that had approached 
or entered the horizon by the time r. The mapping of dark 
properties to evolutionary changes is unambiguous after hori- 
zon entry, when baryonic and photon perturbations are instan- 
taneously impacted by the Newtonian potentials, then reflecting 
the instantaneous overdensity of visible and dark species. The 
remaining, confined to the horizon entry, gravitational impact 
of the dark sectors at a given k in linear theory is also unam- 
biguous (Sec. |2]). 

It is easy to describe the full perturbed linear cosmological 
dynamics (including that of partly polarized photons, baryons, 
realistic neutrinos, quintessence, and other particles or fields) 
by a formalism in which the changes of perturbations in the 



visible sectors are concurrent with the local da rk properties re- 
sponsible for these changes (Bashinsky| |2006l) . The resulting 
description reveals explicitly the objective causal dependencies 
and enables us to map observable features to local dark proper- 
ties more reliably than with traditional formalisms, which lack 
this concurrence. The suggested formalism considers perturba- 
tions of canonical rather than proper distributions. It has addi- 
tional useful technical benefits, for example, note (A) and (B) in 
the caption of Fig. [3] illustrating the corresponding description 
of the acoustic CMB modes. 

Using this formalism, we relate general properties of dark 
perturbations to observable features by tracking the gravitation- 
ally coupled evolution of dark and visible perturbations. The 
advantages of such an evolutionary study over black-box com- 
putation of observables are twofold: It isolates all observable 
signatures of the studied phenomena; it also reveals the mecha- 
nisms that generate these signatures and allows us to judge the 
mechanisms' robustness. 

7.2. Sensitivity of probes 

We categorize the primary probes, responding to the dark dy- 
namics through linear evolution, as either "light" or "matter", 
Sec. [3] Those of the first type (CMB spectra) probe the trajec- 
tories close to null geodesies; of second type (matter transfer 
functions) respond to the metric along the time-like Hubble- 
flow worldlines. 

Primary CMB anisotropies are highly sensitive to the values 
and evolution of the potentials during horizon entry. They are 
only mildly affected by the gravitational potentials on subhori- 
zon scales. (The situation may differ for secondary, nonlinear, 
CMB features.) For example, after recombination the response 
of the CMB angular power spectrum Q to a change of subhori- 
zon potentials is easily quantifiable with eq. ( flTl ), derived in the 
Limber approximation. This response is suppressed relative to 
the contribution (l40l > from the horizon entry by a factor By 
the arguments of Sec. 15.4.31 for linear changes of subhorizon 
potentials this factor is much smaller than unity (~ ri/k). 

The CMB is an excellent probe of the potentials on the hori- 
zon scales at high redshifts. The dark inhomogeneities that en- 
ter the horizon at a redshift z affect most of all the CMB mul- 
tipoles with £ ~ z/20 for the radiation epoch and £ ~ 1\fz for 
the matter era, eq. d42l i. With sufficient angular resolution of 
the detector and reliable subtraction of foregrounds and sec- 
ondaries, the larger number of statistically-independent mul- 
tipoles at higher I's improves the constraints on the parame- 
ters that describe the dynamics at higher redshifts (z, z + Az) as 
yz Az for radiation and \fNz for matter domination, eq. (05]). 

In contrast, matter transfer functions are more sensitive to 
the potential $ at low redshifts. The gravitational impact on 
massive matter at the later times is erased less by the Hubble 
friction. The response of matter overdensity to $(r) at any past 
epoch since the horizon entry can be quantified by a simple 
equation A(6p m / p m )today ~ -tAt k 2 $(r), eq. d24b . The reduc- 
tion of the matter sensitivity to early-time dynamics is mani- 
fested in the prefactor r At. 

7.3. Mapping the dark properties 

The potentially measurable properties of arbitrary dark sec- 
tors may be parameterized by a single function of z for back- 
ground and two (transfer) functions of z and k for scalar per- 
turbations. The effects of dark sectors or modified gravity on 
the metric may be phenomenologically described in terms of 
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{TC, $, \&}. Assuming the Einstein equations, we may instead 
consider either dark species' densities and stresses {p, <5 (c) , a} 
or their dynamical characteristics {w, c\ e , a} (all reviewed in 
Sec.|4]l. The correspondence between the inhomogeneous prop- 
erties of the dark universe and the observable characteristics of 
the CMB and LSS is summarized by Table Q] and the text that 
follows. 

7.3.1. a vs. c e ff 

The measurable dynamical characteristics of scalar dark in- 
homogeneities may be fully quantified by the anisotropic stress 
potential a, eq. ( fT3l . and stiffness ("effective sound speed") cl s , 
Sec. EES 

For adiabatic initial conditions, the gravitational potentials 
$(t,£) and ^(t,Ic) depend on only at the order (Ict ) 2 , i.e., at 
a relatively late stage of horizon entry, Sec. 15. 21 On the contrary, 
anisotropic stress of streaming species affec ts the potentials ear- 
lier, a lready at the order (Jct)°, Sec. 15.11 (Ma & Bertschi nger 
[1991 . This distinction may be given the following intuitive il- 
lustration: The stiffness affects the motion of the dark species, 
hence a certain time is required for the dark species to redis- 
tribute and start sourcing different potentials. On the other 
hand, anisotropic stress is generated without displacing the mat- 
ter. Being itself a source of curvature in the Einstein equations, 
anisotropic stress changes the gravitational potentials earlier. 

There are several observable consequences of the early and 
late influence respectively of a and Cg ff : Both a of freely 
streaming relativistic species and c\ s of a component of dark 
radiation as stiff as quintessence (Cg ff ^ = 1) reduce <j> + \E'. Yet 
the corresponding reductions occur at a different phase of the 
acoustic CMB oscillations. As a result, the CMB oscillations in 
the radiation era are somewhat suppressed by the gravitational 
effect of the anisotropic stress of streaming neutrinos, yet are 
slightly boosted by the gravitationally coupled perturbations of 
tracking quintessence, Fig. fTpl 

Another consequence of the lateness of the gravitational im- 
pact of Cg ff is a considerably larger ratio of the matter response 
over the CMB response to a change of c\ s dark , as compared to 
this ratio for Odark- These results, summarized on the upper half 
of TableQ] should help distinguish an excess of relic relativistic 
particles from a subdominant tracking classical scalar field in 
the radiation era. 

7.3.2. c eff vs. c p 

Two other potentially observable characteristics of the dark 
species that should be distinguished from each other are the 
species' c e ff(fe), considered above, and the velocity c p of the 
wavefront of their localized perturbation. These quantities need 
not be functionally related|3 Together, they are powerful indi- 
cators of the nature of dark sectors. For example, interacting 
relativistic particles whose free-flight time is much smaller than 
the Hubble time have c e ff(£) = c p = 1 /V5; free-streaming rela- 

21 The example of tracking quintessence shows that the reduction of the mag- 
nitude of dark perturbations does not necessarily imply damping of the gravita- 
tionally coupled CMB fluctuations, as sometimes stated. The sign of the impact 
on the CMB anisotropy depends not only on the sign of the change of $ + ^1', 
sourced by the dark perturbations, but also the time of this change. 

22 Since c 2 ff (z,ir) and cr(z,k) provide the most general parameterization of 
the observable dynamical properties of the dark perturbations, these functions 
jointly should encode all the observable signatures of c p . Nevertheless, since c p 
has an important kinetic meaning and maps to a sharp non-degenerate signature 
in the CMB spectra, it is worthwhile to consider this quantity on its own, and 
to compare it with c e g(k). 



tivistic particles have c e ff(k) = 1 / \/3 while c p = 1 ; quintessence 
is characterized by c e ff(£) = c p = 1 . 

Observationally, the phase of the CMB acoustic oscillations 
in both temperature and E'-polarization power spectra or their 
cross-correlation is shifted if and only if c p of any of the dark 
component in the radiation era exceeds the acoustic sound 
speed c s ~ l/y/3, Sec. 1531 dBashinskv & Seliakll2004l) . Any 
such dark component contributes to the phase shift by an easily 
calculable amount. Importantly, the additive shift of the acous- 
tic peaks for adiabatic perturbations is nondegenerate with any 
of the standard cosmological p arameters or with the shap e of 
the primordial power spectrum dBashinskv & Seliakll2004l) . 

Thus for the robust knowledge of the nature and kinetics of 
the dark radiation (encompassing neutrinos, possibly other light 
particles, and early quintessence) both c e g and c p should be tar- 
geted by experimental strategies and data analyses. In particu- 
lar, the determination of c e fr appears most promising from the 
comparison of the CMB and LSS power. On the other hand, 
c p is probed increasingly best by the CMB spectra extended 
toward higher ts, wit h the polarization autocor relation (EE) 
being the most crucial (Bashinsky & Seliak 2004]). 

7.3.3. $anaf$ + 4' 

We can probe the inhomogeneous properties of the invisi- 
ble universe more model-independently by constraining metric 
perturbations directly. Such constraints have a greater validity 
than general relativity, implied for deducing the dark dynami- 
cal properties. They can be placed meaningfully under a weaker 
assumption of only the visible sectors obeying the equivalence 
principle, better constrained for them by terrestrial and solar- 
system probes. 

In this perspective, the discussed probes of the scalar compo- 
nent of anisotropic stress a can be viewed as the direct probes 
of the difference of the Newtonian-gauge potentials $ - 
eq. ( TT~4-b . The primary anisotropies of the CMB are a suffi- 
ciently clean probe of the sum <f> + "J, entirely responsible for 
the gravitational driving of the scalar CMB inhomogeneities 
on all scales and epochs except for a narrow band of redshifts 
around Zdec ~ 1 100. CDM, on the other hand, responds strictly 
to the Newtonian potential $ on all scales and at all times, and 
so does the baryonic matter on large scales after baryon decou- 
pling from the CMB. 

Trivially [eq. (0], the matter growth function is almost di- 
rectly proportional to past $(r) (as previously noted, weighted 
by its conformal time r.) The CMB is also affected by the grav- 
itational potentials. Delay of the decay of <& + "J after the hori- 
zon entry suppresses the CMB power by over an ord er of mag- 
nitud e. [In CDM-dominated limit by 5 2 = 25 times dBashinskv! 
2006) — more than by a factor 2 2 = 4 which appears in the 
apparent description of the Sachs-Wolfe effect in terms of the 
proper Newtonian perturbations.] 

This suppression of the low CMB multipoles is physical: It 
is absent in models where the metric in the matter era is unper- 
turbed. The suppression is diminished in the models in which 
the metric is perturbed less than in the ACDM scenario. The 
prominent suppression of the CMB power at I < 100 by the 
CDM potential is one of the primary reasons that models with- 
out dark matter provide poor fits to the CMB data. The sup- 
pression of the CMB temperature anisotropies for I < 100 must 
not be explained as a "resonant self-gravitational driving" of 
radiation perturbations at / > 200: This effect is caused by and 
probes the inhomogeneities of matter after equality and not the 
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NOTE. — Summary of the discussed properties of the dark sectors, the epochs of their observational relevance, and their effects on the CMB power spectra and on 
large-scale structure. 



inhomogeneities of radiation before equality. 

The suppression of Q for I < 100 severely restricts the 
alternatives to CDM and the models of dark energy which 
reduce metric perturbations at any redshift in the matter 
era. Examples of such mechanisms are contribution of 
quintessence to the density in the matter era, i nteraction or uni 



fication of dark matter and d a rk energy (e.g.JWe tteric 



Perrotta & Baccigalupi 2002; G. R. Farrar, P. J. Peeble 



ch 

es 



1988 



2004 



Catena et al.l l2004; Scherrer 2004), or MOND-inspired alterna- 
tives to dark-matter dMilgrordll983HBekensteinll2004l, . 



minimal physics and therefore should be considered for exper- 
imental constraints whenever possible. Second, GR remains 
falsifiable by other effects; in particular, by the violation of the 
equivalence principle by apparent dark dynamics. This may be 
manifested in the strong dependence of the dark dynamics on 
the visible matter (Sec. 16. U . and in the signatures of effective 
or real superluminal dark flows (e.g., Sec. I6.21 i. GR can also be 
falsi fied by nonstandard phen omenology of gravitational waves 
(e.g. lKahva & Woodardll2007l and Sees. |6\2l and RDl 



7.4. Modified gravity 

Many authors have suggested that modification of gen- 
eral relativity on cosmological scales is the cause of the cos- 
mic acceleration (for rec ent reviews see ICooeland et aDl2006t 
Noii ri & Odintsovl2007l) or even of the apparent m anifestations 
of dark matter jMilgromlll983l;lBekensteinll2004t) . In Sec. we 
consider the phenomenology of typical models of modified 
gravity (MG) that retain the equivalence principle for the vis- 
ible sectors. We show that in these models all gravitational 
impact of the hidden physics can be described within the same 
parameterization schemes of Sec. [4] developed to quantify the 
observable properties of dark sectors that are coupled by gen- 
eral relativity (GR). Indeed, these schemes were restricted only 
by the covariance of the visible dynamics, the assumption of 
the Einstein equations, and the local conservation of the dark 
energy and momentum. However, for any covariant visible 
dynamics, the formal dark energy-momentum tensor d46| i that 
is missing in the Einstein equations is covariantly conserved 
automatically (Sec. |6). Thus all observable signatures of MG 
can be mimicked by effective dark energy and momentum that 
influence the visible species according to the Einstein equations 
and during evolution are conserved locally. 

Particularly, the nonstandard structure growth or ratio 
that are predicted by any MG model of the considered broad 
class can in principle be reproduced without violation of the 
Einstein equations by sufficiently peculiar dark dynamics. Nev- 
ertheless, first, such signatures would still signal some non- 
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APPENDIX 

CMB SENSITIVITY TO (5$ ON SMALL SCALES 

We quantify the CMB sensitivity to the metric inhomo- 
geneities on subhori zon scales and at late times using the Lim- 
ber ap p roximation (iLimberl 1 19541). ap p lied t o the CMB by 
iKaiseri dl984l fl992h and iHu & White! (119961) . The Limber 
approximation for the CMB power spectrum Q assumes that 
the change of the source S in the line-of-sight integral (fJTJ is 
negligible over the wavelength of a typical contributing mode 
k~ l ~ rji. If St is a temporal scale over which the source 
changes by an order of unity then the above condition is equiva- 
lent to k St ^S> 1 . In the Limber limit, C/ is primarily contributed 
by the modes with |k-n|<5T < 1, while positive and negative 
contributions t o AT /T from the p eaks and troughs of the other 
modes cancel dHu & Whitdl 19961) . 

The ISW contribution to the anisotropy source d32l is 

Sisw = + (Al) 
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Mapping the Dark Kinetics 



This direction-independent source in the Limber approximation 
jives the following ad ditive contribution to the power spectrum 



gives tne roiiowing aaaiti 
(lKaiseri[T99i rHu 2000): 



2tt 2 

6Q^— j drrA 2 m {k = £/r), 



(A2) 



where A? c 



1 isw(^) = k Pisw(k) / (2ir ) is the dimensionless power 
spectrum of 5isw(k), evaluated at a time r(r). Thus the ISW 
contribution is 

l7r 1 - ' -t/r\ (A3) 



SO 



dr r g 2 (<£> + 4>) 2 A 2 c (k-. 



where $ and <!/ now denote the potentials' transfer functions, 
i.e., the potentials in the perturbation modes that are normalized 
to a unit primordial curvature perturbation, Cin(k) = 1. 

We now quantify the sensitivity of the power spectrum to a 
change of the potential (5$ on subhorizon scales over a time 
interval St ^S> r/L For an estimate, in eq. (IA3b we set g « 1 at 
the considered late times, neglect the variation of A 2 {£/r) for 
the nearly scale-invariant primordial power, and evaluate the 
remaining integral as 



Then 



SC, 



2n 2 A 2 (£/r) 
I 2 



[8(<S> + ^)Y 



£St' 



(A4) 
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